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HIGH-SPEED SLIP FLOW BETWEEN CONCENTRIC CYLINDERS 


By T. C. Lin and R. E. Street 


SUMMARY 


Schamberg was the first to solve the differential equations of 
slip flow, including the Burnett terms, for concentric circular cylin- 
ders assuming constant coefficients of viscosity and thermal conduc- 
tivity. The problem is solved for variable coefficients of viscosity 
and thermal conductivity in this paper by applying a transformation 
which leads to an iteration method. Starting with the solution for 
constant coefficients, this method enables one to approximate the solu- 
tion for variable coefficients very closely after one or two steps. 
Satisfactory results are shown to follow from Schamberg' s solution by 
using his values of the constant coefficients multiplied by a constant 
factor T| , leading to what are denoted as the effective coefficients 
of viscosity and thermal conductivity. . 


INTRODUCTION 


The fact that a gas is not a continuum but actually a collection 
of molecules in rapid but random motion has begun to have more and more 
importance in the aerodynamics of high-speed flow. This is due to the 
expectation that flow through wind tunnels at low pressure or flight 
at extremely high altitudes will not be amenable to analysis using 
classical fluid dynamics. When the mean free path of the molecules l 
is negligible compared with the macroscopic dimension L, which may 
be wing chord, tunnel diameter, and so forth, the classical picture 
should hold as the molecules are so tightly packed together the gas 
behaves just like a mathematical continuum. Defining the ratio z/l 
as the Knudsen number Kn, it is a measure of the degree .of gas rare- 
faction. In terms of the better known parameters Reynolds number Re 
and Mach number M, the Knudsen number is proportional to M/Re. 

Hence,, telt hough not a new parameter, it is a convenient one to use 
when i 1 degree of rarefaction of the gas is of interest. 

Gas dynamics is simply the ordinary continuous -flow regime or 
Clausius gas regime for which the ordinary Navier-Stokes equations 
together with the condition of no slip on the boundaries are valid, 
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and the Knudsen number is extremely small. If the gas becomes more rare- 
fied and the Knudsen number increases, the effect of slip along the 
boundaries becomes noticeable, although the Navier-Stokes equations - 
remain valid so long as the Mach number remains small. This phenomenon 
has been known for over 75 years and has been the subject of an exten- 
sive 'study by physicists. Tsien (reference l) has summarized this work 
very well. During this same period of time the solution of Boltzmann's 
integral equation by Enskog and Chapman, along lines laid down by 
Hilbert, has led to the distribution function for a nonuniform gas as 
an expansion in powers of the Knudsen number. This approach yields the 
equations of flow in successive orders of approximation, the first 
order being the Navier-Stokes equations, the second order, the Burnett 
equations, and so forth. The third-order approximation has never 
been carried out and the expected complexity of the result does not seem 
to make the attempt worth while, especially as the restrictions on the 
properties of the gas itself are not strictly valid. Chapman and Cowling 
(reference 2) have presented this theory in their well-known treatise. 

Tsien (reference l) presented the Burnett equations of motion and 
pointed out that unless the product Mach number times Knudsen num- 
ber (MKh) was significant any problem in flow could be theoretically 
solved using the Navier-Stokes equations. The question of the proper 
boundary conditions when the higher -order Burnett terms are included was 
raised by'Tsien but not answered until 2 years later when Schamberg, one 
of Tsien' s students, showed in his doctor's thesis (reference 3) that 
the number of boundary conditions required for the Burnett equations is 
the same as for the Navier-Stokes equations.! While being the same in 
number, the Schamberg boundary conditions are considerably more complex, 
being also expansions in powers of the Knudsen number. The first approxi- 
mations for the slip velocity and temperature jump remain essentially the 
same as those used by the physicists in their treatment of low-speed slip 
flows (reference k, ch. 8). The second approximation which is required 
when used in conjunction with the Burnett equations for high-speed flows 
is new and, like the Burnett terms, of considerable complexity. 

The Burnett equations and the Schamberg boundary conditions apply 
to the domain of high-speed slip flow that the aerodynamicist expects 
to enter first when he leaves the domain of classical gas dynamics. 

Their great complexity discourages expectation of a theoretical solu- 
tion of any practically important problem. Hence, the solution of any 
problem, even trivial so far as flows go, is difficult, but if the 
problem can be set up experimentally, the attempt would seem worth while, 
in order to determine the validity of the expansions in powers of. Knudsen 
number and a possible delineation of the dividing line between g^3 
dynamics and slip flow. 


-'-The correctness of Schamberg' s boundary conditions is not universally 
accepted, but they are the only ones. proposed so far. 
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Experimental results have now been obtained for flow past spheres 
(reference 5 ) and transverse flow past circular cylinders (reference 6 ) 

— both- -o-f— which -are-a-lmost— insurmountable— theoret-icaT -problems-? — Schamberg 
( reference 3 ) solved the plane Couette flow problem and the problem of 
the simple rotation between two coaxial cylinders. In order to linearize 
his equations he assumed that the coefficient of viscosity and the coef- 
ficient of heat conduction of the gas were absolute constants. Lin 
(reference 7 ) removed the restriction of constant values of these coef- 
ficients and recalculated the plane Couette flow for a perfect gas with 
constant specific heats and constant Prandtl number but with the coef- 
ficients of viscosity and heat conduction varying 'as a constant power 
of the absolute temperature. The present investigation does the same 
for the flow between coaxial cylinders. In contrast with the plane 
Couette flow problem the flow between two coaxial cylinders rotating 
relative to each other seems of more than academic interest, since an 
experimental check is quite possible and no doubt will be performed in 
the near future. 

j\ 

The problem is set up in its general form assuming only that the l 

Burnett equations and Schamberg boundary conditions are valid and the ' 

flow is steady and stable. Thus the streamlines are circles and only i 

the flow in a single plane normal to the cylindric axis need be con- | 

sidered. Whet her such a flow can be stable at high rotary speeds is \ 
outside the domain of the method us edTEer'e ., A consideration of the ~ 
stability criteria based upon the Navier-Stokes equations with slip at ; 
the boundary is to be found in reference 8 . While special laws of ■[ 
dependence of. the gas properties are assumed, the method is theoreti- 
cally possible for 0 other laws as well as for variable specific heats 
and Prandtl number . 

This investigation was carried out at the University of Washington 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


FUNDAMENTAL EQUATIONS AND EXPRESSIONS FOR STRESS 
TENSOR AND HEAT-FLUX VECTOR 


It is convenient to start from the general equations of the mean 
motion of a fluid all of whose physical properties vary; in Cartesian 
tensor notation these equations are (reference 2 , 7 , or 9 ): 

The continuity equation: 



P 


8u^ 

8x^ 


0 


( 1 ), 
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The momentum equation: 



• dp.. 

pFi + TT-ii = 0 


dx 


J 


The energy equation: 


P 


DE 

Dt 


+ P. 


ij 


duj 8 qi 

8 Xj 8^7 


( 2 ) 


( 3 ) 


where p and E are, respectively, the density and the internal energy 
per unit mass; x A is the Cartesian coordinate in the physical space; 

u i> Fj^, and q ± are, respectively, the components of the velocity of 

the fluid mean motion, the external force per unit mass, and the heat- 
flux vector in the direction; P^ is the component of the pres- 

sure tensor; and 


J. s A * u a 

Dt dt dx^ 


(*) 


is the comoving time derivative or time derivative following the motion 
as in hydrodynamics. (See appendix A for definitions of all symbols.) 
The summation convention, summing over repeated subscripts, is used. 

These general equations can be derived directly from Maxwell' s 
equation of transfer by making use of the properties of the summational 
invariants for molecular encounters without determining the form of the 
molecular-velocity-distribution function. The more convenient Cartesian 
tensor notation is used rather than the vector-dyadic notation preferred 
by Chapman and Cowling (reference 2, pp. 51 and 52). 

In terms of the stress tensor t^, which is defined by 


P-i i = pS-i < + T 


■1J 




ij 


( 5 ) 


(p being the hydrostatic pressure and S-^j, the unit tensor), the 

momentum equation, equation (2), and the energy equation, equation (3), 
take the following forms, respectively: 
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DUi 


St, 


• Sp Ul ii 

p - sr - pF * + S- + srf ■ 0 


( 6 ) 




Su i 

U 8x~ + 

tJ 



= 0 


(7) 


On using the continuity equation, equation (l), and the first law 
of thermodynamics 


dQ = dE + 



( 8 ) 


together with the definitions dS = dQ/T and H = E + p/p, Q, s, and 
H being the heat received, the entropy, and the enthalpy per unit mass 
of the gas, respectively, it is easy to show that 


DE 

P — + 
Dt 


d u i 

c>x . 
1 


= pT 


DS 

Dt 


DH 


Dt 


Dp 

Dt 


(9) 


and energy equation (7) becomes 


or 


_ DS ^ u i 

Dt 8x . 

J 



pSs.auT,, 

Dt Dt ^ 8x^ 8x^ 


0 


(10) 


( 11 ) 


Adding to equation (ll) the product of u^ and equation (6) and 
considering that Tjj is a symmetrical tensor, another form of the 
energy equation is obtained 




P u i F i 



( 12 ) 
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Equations ( 3 ) > ( 7 )> (10), and (ll) are the most familiar forms of the 

energy equation. Neglecting terms containing the external force 

and setting H equal to c p T, equations (6) and (12) yield the general 

momentum and energy equations given by Tsien (reference 1) and used by. 
Schamberg (reference 3). The sign of T.y, being the same as that used 

by Tsien and Schamberg* is opposite to the usual convention. 

It is only through the expressions for the stress tensor and 

the heat-flux vector q^ that the above momentum and energy equations 

depend on the form of the molecular-velocity-distribution function. 

Let r T^ and j> q i denote the rth-order approximations to the stress 

tensor and the heat-flux vector q i , respectively, and write 

r T ij * )L T iJ <n) (13) 

n=0 

and 

=11 9i (n) (l 1 *) 

.n=0 


where 


t. .( n ) 
ij 


and q. 


(n) 


are the nth-order corrections to 


0 T i j 


and _q., respectively. Then the first-order approximation to the 

mole'cular-velocity-distribution function, that is, the Maxwellian dis- 
tribution, gives (reference 2, pp. 112, 122, and 123) 


0 T ij 


= T.. (0) - 0 

1J 


( 15 ) 


and 



(°) = 0 


(16) 
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which, together with equations (6) and (7), yield the Eulerian equation 
of motion 


P 



- pF ± + $ 2 - = 0 


(17) 


and the corresponding energy equation 


p de + 

Dt 


8 u^ 

8x^ 


= 0 


(18) 


The second-order approximation to the molecular-velocity-distribution 
function gives (reference 2, pp. 112, 122, and 123) . 


l T iJ ' T iJ 


( 1 ) _ 


= -2p 


^i 

Bx j 


(19) 


and 




(20) 


where p and X. are the coefficients of viscosity and thermal conduc- 
tivity, respectively, and is the nondivergent symmetrical tensor 

J 

du-j 

associated with the tensor - — that is. 


^ u i _ l/^ u i + duA 1 

2 U Xj 8 Xi y 3 a Xk °iJ 


5 . 


In general, any tensor A-y with a bar over it has the following 
meaning (reference l) 


A ij “ 2 ( A ij + A ji) 3 


A kk^ij 


( 21 ) 



8 


NACA TN 2895 


Substituting equations (19) and (20) into equations (6) and (7) 
yields the conventional NavieivStokes equations 


Du-: A-r, £ / 8u,- 1 

c 5T - pF i + SI - 2 ^ s£) ' 0 


dx A Sx j> 


( 22 ) 


and the corresponding energy equation of viscous flow 

chi.: du,- du.* 


de A du i a du i du i a A a r \ _ n 

Dt p a Xi a X j a X j a Xi y a Xi j 


(23) 


From the definition of the nondivergent symmetrical tensor, equation (21), 

M S) ' 3 si( M sj) * 2 ir } (( sj) 


and 


a Ui a Ui a Ui au^ 
a^r a^r “ a^r a Xj 




au 2 ^ 2 


L SXq 




au 3 \ 2 


du- 


a Ul ' 2 


, a Xo ax, 


. a X o a Xn 


fou-L au 3 \ 2 /au 2 auA 2 / au 3 

^a X3 a Xl y Va x ’i ax 2 y . ya X2 


au 2 x2 
ax, 


^ 0 


Making use of these relations, it is readily seen that equations (22) 
and ( 23 ) check with the momentum and energy equations of viscous flow 
(reference 10). In the present notation the dissipation function is 
simply 


au^ a Ui au-^ au^ 

* " ^3~ 3x j 


( 24 )- 


which is always positive and unaffected in form by the fact that the 
coefficient of viscosity p is a variable. 
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The third-order correction to the molecular-velocity-distribution 
function, _as given by Burnett (reference 2 or 11), yields °the second- 
order corrections to the stress - tensor - arid the~heat-f lux vector - for both ~ 
spherical and Maxwell molecules. These corrections are accurate to terms 
of order (p/p) 2 . In Cartesian tensor notation they are (references 1, 

3, and 7) 2 



K-, 


p.2 du k 8 Ui 
P Sx k ^ x j 


k 2 h! 

mi 


_L/f. 

8x ± l j. p Sxj j 8xi 8x k 8x k 8x, 


K 3 + K, ttL |2L + 

J pT ox£ oxj ^ ppT Sx^ Bxj 


p 2 8t 6t + K p 2 ^ u i ^ u k 

0rj, 2 dx^ dxj 6 p 8x k 8xj 


(25) 


q (2) . ai |s_ 

1 ■ L pT ox^ oxi 


e 2 


£ 

2 

8 

4 ^ 

.2^1*- 

PT 

3 

8 Xi 

\ 5x j/ 

8 Xi 8 Xj 


(' 


e 3 <*_ . e 4 e! a * 95 e! | r\^i 

PP 8xj p dxj ? pT 8xj/8x^ 


( 26 ) 


^The last term within the brackets in equation (25) differs from 
the one given in the references quoted. The correct form for the Burnett 
terms has been given in reference 9. This error was pointed out to the 
authors by Prof. S. A. Schaaf arid confirmed by Prof. C. A. Truesdell and 
Mrs. C. S. Wang Chang. 
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The K' s and 0's are pure constants. Their correct values are given 
below (references 7 and 9): 


Maxwell molecules 


Rigid elastic spherical molecules 


v _ h(l T dp\ 

K 1 - 3V2 ' ir dT/- 

K 2 = 2 

K3 = 3 

K4 = 0 

v _ o T du 

*5 ■ 3 iT dT 

K 6 = 8 


>(27) 


01 

0 2 

0 3 

04 


15 

4 


(1 1 

\2 " p dT/ 


-i£ 

8 


= -3 


= 3 


M29) 


= M + £ ^ 

\4 H dT/ 


K i ■ 1(5 - ? ®) x 1 - 011 ‘ 


K 2 = 2 x 1.014 


K 3 = 3 x O.806 


Kj^ = 0.681 

*5 = 3 ~ dT x °' 8 ° 6 _ °*990 

i 

Kg = 8 X 0.928 

e, = 15(1 - 1 St) x 1.035 

1 4 \2 n dT/ 


0o = 


0O = 


3 

04 


^ X 1.035 


-3 x 1.030 
3 x 0.806 


(28) 


>(30) 


= 3(^2 x 0.918 + I dix x 0.806^ 

\4 h dT / 


0.150 


It is noted that the values of 0 2 and 0^ for Maxwellian molecules 

given above are different from those given by Chapman (reference 2, 
pp. 267 to 270). These corrections together with the values of 0 for 
rigid elastic spherical molecules are due to. Wang Chang and Uhlenbeck 
(reference 9). For ordinary gases (T/u)(dp/dT) has a value lying 
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between l/2 and 1. It follows from the corrected expression for 9 
given above that all the coefficients and are less than 117/4 

instead of 45/4 as given by Chapman and Cowling (reference 2, p. 270). 


With the expressions for the stress tensor and the heat-flux vector 
accurate to the second order, the momentum and energy equations, equa- 
tions (6) and (7), become 



pFi 


+ 




(3D 


Dt dx.j. dxj dxj dxj l chew 1J 8xj d Xi 


0 ( 32 ) 


where t^ 2 ) and q^ 2 ) are given by equations ( 25 ) and ( 26 ), respec- 
tively. Equations (3l) and (32) are the momentum and energy equations 
for slip flow and they reduce to the Navier-Stokes equations and the 

corresponding energy equation on neglecting terms containing 
and q^ 2 ). J 


From equations (19), (20), (25), and (26) it is seen that the ratio 
of a typical term of t^ 2 ) to or qj 2 ) to ± q ± has the same 

order of magnitude as either |iU/pL or |a 2 U/X.pTL, where U and L are 
the characteristic velocity and length of the flow, respectively. Making 
use of the relations 



and 


a 2 = 7RT 


(33) 


(34) 


where R is the gas constant, Cp, the specific heat at constant pres- 
sure, 7, the ratio of the specific heats, and a, the adiabatic speed 
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of sound, together with the definitions that the Reynolds number 
Re = pUL/p, Mach number M = U/a, and the Prandtl number Pr = Cpp/k, 
it follows that x 

— U = (7 - l)Pr M 2 /Re (' 

kpTL 

From equation (34) and the' perfect-gas law v 


it also follows that 


p = RpT 


— = 7M 2 /R€ 


From the kinetic theory of gases (reference 4, pp. 50 and 147) 

p/p = a 2 /7 = jic 2 /8 


n = 0.499 pcZ (3S 

where c" is the mean molecular speed and Z is the mean free path of 
the gas. From equations ( 38 ) and (39) 

! = ^598- r , - 2 M/ap <1,c 

Hence the Knudsen number is 

Kn = Z/L = 1.2567 1 / 2 M/Re (4l 

Substituting equation (4l) into equations (37) and (35)> 

— = 0.7967 1 / 2 M Kn (42 

pL 


llfiL.o 


• 796 (/7 Tf* M 
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For ordinary gases the Prandtl number is approximately unity and 
~ 1 ^ therefore pU/pL and p^u/lpTL have the same order of mag- 
nitude as M Kn. Since the rarefaction of the gas increases with the 
Knudsen number, it is evident that for the high-speed flow of a rare- 
fied gas the second-order terms r^ 2 ) and q.^ 2 ) Q f stresses 

and the heat flux become relatively important. For 7 = 1.400 and 
Pr = 0.750, it follows from equations (4l), (42), and (43) that 

Kn = 1.486M/Re, pU/pL = 0.940M Kn, and p 2 U/XpTL = 0.201M Kn. 

According to Burnett's expression for the molecular-velocity- 

(3) 

distribution function, Tjj , the third-order corrections to the 

stress tensor, will contain terms of the form (p/p)3^du the 
ratio of which to 2 . T ij ^ as same order of magnitude as 



M 2 Kn 2 


(44) 


Similar terms apply for 'IV'* Hence the slip-flow equations, 

equations (31) and (32), cease to be valid if, for a given Mach number, 

the gas is so rarefied that M 2 Kn 2 is not negligible compared with 
unity. 

4 . 

The particular problem to be considered in this investigation is 
the slip flow between concentric cylinders. This problem has recently 
been solved by Schamberg (reference 3, ch. VII) for the case of constant 
coefficients of viscosity and thermal conductivity. The present inves- 
tigation extends Schamberg' s solution to include the effect of variable 
coefficients of viscosity and thermal conductivity. 

Assume that the rarefied gas is confined between two concentric 
cylinders. The inner cylinder, having radius a and the uniform tem- 
perature T wa , is rotating at constant angular velocity cn^, its sur- 
face velocity being denoted by U = aai tfa ; whereas the outer cylinder, 
having radius b and the uniform temperature T v ^, is held fixed in 

space. The flow field is conveniently described by the cylindrical 
polar coordinates r, 0 , and z, with the z-axis as the axis of the 
cylinders. As sum ing that the flow is two-dimensional and steady and , , 
that the external force~caiT.be ’ neglected, • ‘ 
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u 


z 


= 0, 




Fi = 0 


It follows from the symmetry of the problem that 


( 45 ) 


U r = °> h = 0 (46) 

00 

where u z and Uj. are the velocity components in the z and r direc- 
tions , respectively. 

The appropriate equations of motion are obtained by expressing the 
continuity, momentum, and energy equations, equations (l), (31), and (32), 
in plane polar coordinates. This is easily done by making use of the 
formulas in general orthogonal coordinates (references 12 and 13). In 
view of equations (Vpj and (h6) it is found that the equation of con- 
tinuity is automatically satisfied. The momentum equations in the tan- 
gential and radial directions are, respectively (see appendix B), 



The energy equation becomes 



In the above equations t r j is the viscous shearing stress, T rr and 
T 00 are the normal stresses in the radial and tangential directions, 
respectively, and q r is the radial component of the heat-flux vector. 
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The explicit expressions for the required components of the viscous 
stress tensor and the heat-flux vector are obtained by transforming the 
general _expressions, equations _(23.)- and_.(26),_ respectively, into plane 
polar coordinates and making the reductions required by equations (45) 
and (46). This is done in appendix B. The results are given as follows. 



(50) 


rr 






3 /\ T / 3 4 pT dr dr 3 rp dr 




i KoR 1 ^1 
3 r drj 


(51) 



3 3 dr 2 



2 K « ® T 

3 3 r drj 


(52) 


<3r = 



(53) 


The expressions for T rr and T00 given in equations (5l) and (52) are 
slightly different from those used by Schamberg (reference 3, PP. 152-153) • 
This is due to the correction of the last term in the brackets of equa- 
tion (25) for the Burnett terms in the stress tensor. 
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Substituting equations (50) , (51 ) 7 (52), and (53) into equations (47), 
(48), and (49) and using equation (27) for the values of the K' s for 
Maxwell molecules give 



T 1 

where p = — Thus one has three equations, equations (54), (55), 

and (56), which together with the perfect-gas law, equation .( 36) , are 
used to find the four dependent variables u^, T, p, and p as func- 
tions of r. 


BOUNDARY CONDITIONS 


It is seen from the previous section that the introduction of the 
higher -order approximations to the stresses and the heat flux results 
in an increase in the order of the momentum and the energy equations 
of the fluid mean, motion. For instance, the first-order approximations 
to the stresses i T ij, equation ( 19), lead to the conventional Navier- 

Stokes equations, equations (22), of viscous flow, which are partial 
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differential equations of the second order; while the second-order 
approximations to the stresses lead to the third-order partial differ- 
ential equations, equation (3l)> with the expression for T given 

by equation (25). .Since.bhe relative importance of the higher-order 
terms of the stresses and the heat flux increases with the rarefaction , 
*Pf~ttie 'gas, this leads one to the expectation that the number of bound-. ' 
ary conditions required for the complete evaluation of a slip-flow j 

"problem should likewise depend on the degree of. the rarefaction of the J 
"gas, however, it was shown by Schamberg (reference 3) on both physical"]; 
ana ' mat hemat i c a 1 grounds that the number of physical boundary conditions j 
required for a slip-flow problem is effectively the same as that for 
the corresponding flow in the realm of gas dynamics. ^ 

In slip flow, as in gas dynamics, the condition of zero relative 
normal velocity at the boundary still holds, but the relative tangential 
velocity at the boundary is no longer zero and the gas temperature dif- 
fers from the wall temperature.. These are known as the "slip velocity" 
and the "temperature jump," respectively. 


The expressions for the slip velocity and the temperature jump at 
low Mach number were investigated by Maxwell, Millikan, Smoluchowski, 
Khudsen, and others (references 1 and 4). If x and z are the dis- 
tances tangential and normal to the wall, respectively, u and u^. 

are, respectively, the velocity of the gas and of the wall in the 
x direction, and T and T w are, respectively, the temperature of 

the gas and of the wall, one has from the kinetic theory of gases (refer- 
ence 4) the slip velocity 


and the 


< u) z=0 - “w * 
temperature jump 




3jj_/8t\ 

4 pT\8x/ z=0 



(57) 


(58) 


where a is Maxwell’s reflection coefficient, c v , the specific heat at 
constant volume, a, the accommodation coefficient, and 2, the mean free 
path. 

For slip flow at high Mach number the above expressions give no 
longer a true description of the physical relations and higher-order 
approximations to the slip velocity and the temperature jump must be 
used. A general method for the calculation of approximate expressions 
for the slip velocity and the temperature jump, to an arbitrary degree 
of approximation, is given by Schamberg (reference 3). The method 
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applies the laws of conservation of mass, momentum, and energy to the 
infinitesimal layer of gas adjacent to the solid surface, referred to as 
the "sublayer," and uses the nonuniform molecular -velocity-distribution 
function. 

The first approximations for the slip velocity and the temperature 
jump thus obtained by the use of the first approximation to the velocity 
distribution for a monatomic gas of Maxwell molecules agree with the 
results given by previous investigators for low-speed slip flows. 


The second approximations to the slip velocity and the temperature 
jump, which are required in conjunction with the second approximations 
to the viscous stresses and the heat flux for high-speed slip flow, are 
given as follows (reference 3): The slip velocity 





9 r _8t dw 8_ RT dp/d u j + aw_\ _ 

8 8 z 8 xj 15 P dxj/ 

£ lj)] + (^ <RI)1/2 0 * Ij) ♦ 

i M sj - 

a 1 (RT) 1 / 2 fci + /log H.\ _ 

l8z dxjyDt \ e P / 


1 R D /^T \ . RDT 8T_ 9 R Dp 8 t_ 

8 2 T Dt dxj 8 P Dt 8xj 

4 R DT dp_ $p u j + dw \dw 

5 p Dt dxj 6\dz ^Xjlpt 


( 59 ) 
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where j 
T(x,y,z) 


=1 a nd 2, and the temperature jump 


z — >0 


T » - + (;) 2 [ lT p + 1 ) 

I T I 


J | 8w 


e 3 T(ET)1/2 ( 2 ^ + sflj * 
e * (BT)1/2 l^r * t) * |j 

i*"^l^ + sWs^- 


J J 


e p 8 t_ 8t_ + e r (8t\ 2 , 1 RT 8t 8p f 

' cbcj ® \^ z / 8 p dxj 8xj 

e ^p_ <*L _ _L RT . _ J_ RT ^ c>hl 

^ p 8x j 3x j l4 ^ z 2 1^ p bz dz 

/fci\ 2 J. 1/DT\ 2 . 1 DT D L nn p\ a 
(p) j 9 l\pt) 2 Dt Dty 1 g e pj + 


'10 


( RT \l/2 1 DT 8 t / rtV 1/2 D /8t^ 
(RT) T Dt" "Sz 6 5^ RT) DtteJ 




3 e 3 


1 up Dw l8 Dw 
(RT )l/2 Dt &t " T 87 Dt 


where j = 1 and 2. 
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In equation ( 60 ) but not in equation ( 59 ) the summation convention 
over two indices is used. This notation has the meaning: 


X x = x > x 2 = y 


u-^ = u ; U2 = v 


(61) 


so that 





dh/c)u + 
8x \^z 





(62) 


The time derivative D/Dt, when expressed in the above notation, is 


— = — + u ^ 

Dt 8t W J 8xj 


( 63 ) 


All of the derivatives in equations ( 59 ) and ( 60 ) are to be evaluated 
at a. point (x,y,z) as z — >- 0 . These boundary conditions are appli- 
cable provided that the. pressure level and the motion of the gas are 
such that the second approximations to the stress tensor 2 T ij and 
the heat-flux vector 2^1 are applicable. This means that the rela- 
tion M^Kn^ « 1 holds. 

The values of a^, the b's, c^, and the e' s in equations ( 59 ) 
and (60) are given on the following page (reference 3):^ 

Alchamberg' s values for these constants have not been changed to 
agree with the corrected values of the K' s and 6 ' s given by equa- 
tions (27) and (29) because the corrections are negligible to the order 
of approximation used later on. 
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In general c-^ is used instead of c-^. According to the kinetic 
theory for monatomic gases (reference 2) .7 = 5/3; X/^ c v = 5/ 2 > and 

c^ f reduces to c^. The values of a for air vary from 0.79 to 1.00 
while a for air lies between 0.88 and 0.97 (reference.!). 


To obtain the boundary conditions for the concentric -cylinder flow, 
it is first necessary to express the general boundary conditions, as 
given by equations (59) and (60) in the Cartesian coordinate system, . in 
terms of the polar coordinate system. Following Schamberg (reference 3) 
closely, let x and y be the Cartesian axes associated with the polar 
coordinates r and 0 and the auxiliary coordinate systems x a , z a 

and x^, by the equations 


x = a + z a = r cos 0 




y = x a = r.sin 0 



( 65 ) 


x = b - z-^ = r cos 0 
y = -x b = r sin 0 




( 66 ) 


The velocity components u a , w a and w-^ of the auxiliary 

coordinate systems are related to the tangential velocity up as shown 
by figure 1 and equations (67) and (68): 


u a = ^ cos 0 
w a = " u 0 sin 0 


( 67 ) 


u b = " u 0 cos 0 
w b = u 0 sin 0 


( 68 ) 
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The partial derivatives vith respect to the auxiliary coordinates are 
expressed in terms of the partial derivatives with respect to r and 

_0 by means of eq ua tion s (69) and (70) » which a re easily obtained from 

equations ( 65 ) and (66), respectively: 


a 

dx a 


sin 0 — + cos 0 — — 
dr r 80 


> 



(69) 


a 

8x b 


-sin 0 — 
dr 


cos 0 — — 
r 80 


— = -cos 0 — + sin 0 — — 
8z b 8r r 80 


(70) 


All of the first- and second-order partial derivatives appearing 
in equations (59) and ( 60 ) for the- boundary conditions can now be 
transformed into polar coordinates by means of equations ( 67 ) to (70). 
After all of the differentiations with respect to 0 have been per- 
formed, 0 is put equal to zero in the resultant expressions, in 
accordance with figure 1. Because of the condition of axial symmetry, 
equation (46), the partial 

respect to 0 all vanish. 

The boundary- condition derivatives for both the convex (r = a) 
and the concave (r = b) cylindrical surfaces are given in the table 
on the following page. The transformations for the derivative D/Dt 
are obtained from equations (63), (45), and (46), using the fact that 
the velocity of the wall u w has the value U at r = a and zero 

at r = b, respectively. 


derivatives of 




T, p, and p with 
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With these substitutions one obtains from equations (59) and (60) 
the following four boundary conditions: 

(l) At r.= a. 




(71) 


(72) 


(73) 


(4) At r = b, 



(74) 


The subscripts a and b denote the evaluation of a particular quan- 
tity at r = a and r = b, respectively. The quantities X a , X b , 

Z a , and Z^, are defined by the following equations: 
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The values of the numerical constants a-^, c-^, e-^, eg, and eg are 

given by equations (64). 

SOLUTION OF CONCENTRIC-CYLINDER FLOW 

From the preceding sections it is seen that the problem of slip . 
flow between concentric cylinders is reduced to that of solving the 
three differential equations, equations (54), (55) > and (56), and satis- 
fying the four boundary conditions given by equations (71) to (74). For 
this purpose it is desirable to introduce dimensionless constants and 
variables as follows: 


k = a/b 


(79) 
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Mua 


Re. 


wa 


Pr. 


wa 


j^va_ 

p a Uh 

Hwa 

^ wa - C Pwa 
\t a 


(80) 


!•* = r/a 

P* = P/P a 

\* = V*™ 



u* = yju 

P* = p/p a 

c p* = c p/ c p 

*1 ^va 

> 

(81) 

T* = T/T wa 

P* = p/p wa 

Pr* = Pr/Pr va 




where 

h = b - a . 3 (82) 

and c p wa > Pwa' and \ia. are the P r °P e:r “ties of the gas based upon the 

wall temperature T wa . Making use of the perfect-gas law, equation (36), 

and of equations (79) to (82), equations (54), (55), (56), and (36) 
become, respectively, 



(83) 




= 0 


(84) 
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d log e P* 
dr* 


1 dp* 
p* dr* 


= 


2 (u *) 2 (1 - k) 2 7 M ^ a 2 

" ' k 2( Ta *)2 Rewa 2 


(n*)‘ 


r*(p *) 2 




o „ 2/u*\ 2 , dV , 1 dT* 3p/dT*\ 2 

27M wa ( r *) 3 ( dr *)2 3 r* dr* T*W*/ 


( a io ge P *\ 

m* d log„ p* 

1 0 1 c V 

\ dr* / 

r* dr* 1 

*— 


dr* 


(1 - k) 2 7M wa 2 d 


i(n*V 


r 2 2/du*\ 2 

[_ 3 7 ^ a \dr*/ 


k 2 p*(T a *) 2 Re W a 2 dr * L P 

4 m 2 /V\ 2 p d^* 1 dT* { 2p/dT*V 

3 7 ^ a \r */ . (dr*) 2 ’ r * ^ T* \dr*/ 


,*\2 


“ “^*( T 
3 dr*\ 


L» a l0 «e p *) 

2 T* d 

log e P* 


\ dr* ] 

3 r* 

dr* 



(85) 


p*T a * = P*T* 

Velocity and Temperature Distribution 
Integration of equation ( 83 ) gives 


( 86 ) 


u* - Ar* - 2Br* 


I 


dr* 


(r*) 3 p.* 


(87) 
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where A and 2 B are constants of integration. For u = Constant, 
U* = 1 and equation (87) yields 


u* = Ar* + - 2 - 

r* 


(87a) 


which is Schamherg' s solution. 

From equations ( 84 ) and (87) 


<r*> V £;(**>■* ^) * Mr - l)Pr va M ua % 2 = 0 (88) 


It is convenient to transform the independent variable r* into £ 
or 5. The latter are defined, respectively, as 


log e r* 1 

t = — = - log r* m = log l/k (89) 

log e l/k m e e 


and 


?-n /*a (9o) 

J o V* 

where 



Equations. (89) to ( 91 ) and ( 79 ) insure that 


1=^=0 at r* = 1 or r = a 


£ = i = 1 at r* = l/k or r = b 
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In terms of | equation ( 88 ) becomes 


P* — ) + Mr - l)Pr wa - 0 (92) 

while in terms of the independent variable £ it is 


if 



Mr - 1) 


Pr. 


wa 




0 




(93) 


Equation (93.) together with equations (90) and (91) gives A(o) = A(l) = 0 
In general A(£) « 1. 

Integrating once 

£ jP - - - ft.. (*?)£*<« d5 

and. again 


T* = D + mC d£ - 1} Pr wa M^ 2 B 2 J £ d tj (k 2 )^ +A( ^ d£ 

(9M 

For constant (i and X, |i* = X* = 1 , q = 1 , £ = and equation (9*0 

reduces to 


T* = D + C log e r* - (7 - l)Pr wa M wa 2 B 2 (r*) 2 
which is Schamberg f s solution. 
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In terms of equation ( 87 ) becomes 


* u * A 

= — = A 
r* 


_ 2Bm J (k2) 


5+a(S) 


d£ 


(95) 


Now let 


' n=l n * j 


(96) 


This series converges for all finite values of m A(£). Making use of 
equation ( 96 ), equation (95) gives 


where 



(97) 


00 

f(0 =H 


n=l 


( -2m) n+1 
ni 


f k 2 ^(A(cTl n d£ • 


(98) 


Similarly , from equations ( 96 ) and ( 9*0 , assuming Pr/c p = Constant, 

T* = D + mC£ + 7 —^ Pr wa M wa 2 B 2 [k 2 ^ + g(£)] (99) 


where 


g(5) = 2m J f U) dg 


(100) 


i = 


Now at r = a, £ = 0, <d* = co a *, and T* = T & *; and, at r = b, 

1, cd* = o>^*, and T* = T^*. Therefore, from equations (97) and (99) 


ay 


- 


= 0) a * - (m a * - cDb*) 


1 - k 2 ^ - F(£) 
1 - k 2 - F(l) 


( 101 ) 
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where 


F(C) = f(0 - f(o) 


( 102 ) 


and 


T* = T a * - (T a * - T h *)t + 

. t ; k 2 ? . k 2? 0U 3 (103) 

[1 - k 2 - Fdl] 

where 

’V 

G(g) = g(U - g(0) - (g(l) - g(oT|c (104) 

It is noted that F(0) = G(0) = G(l) = 0. 

Equations (101) and (103) represent the exact solutions of differ- 
ential equations (83) and (84). Since the latter did not contain any 
Burnett terms, it' was relatively simple to carry out the integration. 
However, A(£) depends on a knowledge of rj and which are func- 

tions of T* and so any numerical solution will reduce to an iteration 
process of approximating A(£) and hence F(£) and G(£)» This is 
also the method of determining the pressure distribution p* from 
equation (85). Before proceeding to p*, it is advantageous to rewrite 
the boundary conditions of the section "Boundary Conditions" in a similar 
dimensionless form since the constants o) a *, 0^*, T a *, T b * i n 

equations (101) and (103) are precisely these boundary values. The 
actual algebraic reduction is carried out in appendix C, leading to the 
following results: 



,1.592a 


V - <°b* 1 - k 

(Ta*) 1 / 2 . k [T - k 2 - F(lT] 




«*e + \ 



(105) 


= 1.592a 1 


K* - %*)( T b*) l/2 k 2 (l - k) 


T a *V 


1 - k 2 - F(l) 


Kh e + kX b Kn e " 


(106) 
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T* = 1 - 0.796c 


1 - k 


-km^Ta *) 1 / 2 


T * - Tt_* + 
a b 


(7 - l)Pr wa - <%*) 


2 1 - k 2 - 2m - G'(0) 
Q. - k2 - F(l)] 2 . 


f Kn e + Za Kn e ^ 

(107) 


T b * =• T wb * + 0.79 6 C;l 


(1 - k ^ Tfr*) 1 / 2 


Ta* - T b* + 


(7 - l)Pr wa “ra 2 ^* " ®b*) 


2 1 - k 2 - 2k 2 m - G* (l) 
[l - k 2 - F(l)] 2 


where 


fKh e + Z b Kn^ 
(108) 


^e = ^a 


(109) 


'W (V 

and X a , X b , Z a , and Zjj are the complicated expressions given in 
equation ( C12) . Actually, since X & . . . \ are multiplied by Kn e 2 


above, it is only necessary to use the zero-order approximation to these 
expressions given by equations (C32), (C33), and (C34). The reason for 
this, of course, is that the differential equations and the boundary 
conditions have been derived only to terms in the square of the Knudsen 
number. Thus, it is probably clearer to write equations (105) to (108) 
in the form (equations (C27)) 

“a* - o 

+ ,u>*Kn + 0 (x> *Kn 2 ) 
la e 2a e / 

“b* = tfV + 

1%* ^e + 2%* Kn e 2 

T a* = oV^ 

+ !T a *Kh e + 2 T a *Kn e 2 ) 

T b* = o T b*(i 

+ iVia>e + aV^e 2 ) 
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which explicitly indicates an expansion in powers of Kh e . The coef- 
ficients h co b*4 h^a*' h^b* (^ = 0, 1, 2) are written out 

in equations (C29) to (C3l) of appendix C. The subscript written in 
front of a symbol thus denotes the order of the approximation in the 
expansion. 


Pressure Distribution 


From equation ( 85 ), neglecting the terms containing 


Vfte„a/ 


which 


P P 

are the order of co or Kn , the zero-order approximation to p* is 



(no) 


The superscript in front of a symbol denotes the order of the approxi- 
mation to the solution of a differential equation. 

From equations ( 89 ) and (90) 


dr* 

r* 



Therefore 


Op* = exp 






( 111 ) 


or 


Op* = exp 


ZS M 2 H*(<o *) 2 

_n Jo T * (k 2)^A(c) 


(112) 


Since equation ( 85 ) does not explicitly contain terms of order Kn, the 
zero- and the first-order approximation for p* can be obtained from 
equations (110), (ill), or (112). From equations (112) and (C29) 
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From equations (85), (C 5 ), and (110 ), log e 2p* = log e °p* + f Kd^ 2 
or • 

V * (1 + 0* If (114) 

with exp Kn^ 2 ) .% 1 + Kh wa 2 where 


0* = 


0 . 6336(1 - k ) 2 

k2 



l(n*) 2 dr* 


r* 0 (p*)‘ 


-s27M wa 


/u* 

n\ r * 


££ W3 — + 

te*J o\ r* 


du*\ 

dr*/ 


+ 



0.6336(1- k) 2 
k2 



dr* 

OP* 


d 

dr* 



2 2 


/u* 

oV r * 


du*\ / u* + du*\ 
dr*/ r* dr*/ 


+ 



is obtained from the second-order terms of equation (85) upon substitu- 
tion of the zero-order approximation therein. 
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In terms of equation (115) becomes 


0* 


0.6336(1 - k) 2 




0^ 


( do)* 4m 


Wo 


0 







d£ + 


2. P J_± J_[i 2 n M 2 (-££ - 6 2 u<J\ 

m Jo O 15 * d £ 0 V *\l m 0 Va <W / oV d ^ n ' 


a 2l (k2) C*A(C) /dM . 3 , (k2) 

m o\ d r / 0 


5+A(5) / dT* ' 

oV d £ > 


d £ 


(116) 


Integrating by parts, using equations (C33) and (C34), and neg- 
lecting A(£), equation (ll6) becomes 


»+ = - 


0.6336 

fo 1 ! 

0^* / 

k 2 (l + k) 2 

iT 

L (p *> a \ 


- 2 S «- ^k^k 2 * 


0^ 


5 


*f 


»— 

d 

1 1 

^ d^ 

b(P*) 2 J 


^ 3N l " 6N 2 k ^ ' 8 n 3- “ k ^j k2 ^ d C + 


6 P ° ( ^ -L|n 1 - 2N 2 k 2 ^' - 2N 1+ k 2 ^k 25 d^ 




(117) 
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where 


N, 


* (1 - T w b*)(i - k 2 )' 2 * {y - l)Pr wa My a 2 {l - k 2 ) 

I ^a 2 ” 


Nk 


, n 3 


= Mwa 2 ®^ - DPr wa - ^ 

= [3 7 + (7 - 

^Mwa 2 ” 


- 7 + (7 - l)Pr wa 


(118) 


As a first approximation put = qT] and assume 


with 


— = Aq + A x k 2? + A 2 k^ 

o(p *r 


Aa = 1 + (2k 2 - 1) 


A 1 = -2k' 


1 ~ o(V) 

(1 - k2) 2 

1 - oM~ 2 


-2 


a 2 = 


(1 - k2) 2 

1 ~ o(Pb *)~ 2 

(1 - k 2 ) 2 


(119) 


(120) 


Equations (119) and (120) insure that qP* = 1 at' £ = 0 and that 
qP* = qP^* and d^p*j d £ =0 at t, = 1. Then from equation (117) 

°* (£) " ‘ AjA 3 + V^ 1 - kH ) * 

( A 2 A 3 + a 1 a 5)( 1 " k 6 ^) + AgA^l - k 8 ^)j 


( 121 ) 
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where 


■ <7 - «»a 2 ♦ (l - T ub *) 


,(1 - k 2 ) 2 


•| ^«»a 2 


A4 ■ (§ r - 8(7 - Dpt^m^ 
A 5 = It r ‘ T (r ' 1)Pr ™] M w= 

a 6 - i ) pr v3 M " a2 


( 122 ) 


Then for £ = 1 




E ■ o(^*) 2 ] B (3 * 


26k 2 - 17k 7 + 


7 - 1 


Pr wa (3 + 2k 2 - ^ k ) + 


7 - 1 ^ 

2 ^"wa 


(1 - k 2 )( 


)(2 + k 2 ) 

m 


0.6336(1 - k) 2 (2 + k 2 ) 


)(V) 


(123) 


The integral in equation (ll4) remains to be evaluated. It is 
necessary to assume some relation between n* and T* which, to avoid 
too much complexity, will be assumed to be the simple relation n* = T*P 
where P is a constant whose value lies between 0.5 and 1.0 (reference k, 
p. 150). Furthermore, in many cases p s: 0.9 (reference l4) so the 
assumption P = 1 or ^* = t* will be made in order to simplify the 
evaluation of this integral as well as to enable a more direct expansion 
of all distributions later on. 
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Equation (ll4), with the aid of equations (101) and (96), then gives 
for 0 = 1 ~ 


2 p* = J7 + 0 t(5)Kn e ^j exp 




2q[l - k 2 - F(lT]' 


jK* - cu b*) 2 ( 1 - k 2 ^) “ 


4m(o> a * - V )|kV + F ( 1 )“ a * _ + 


\*\* + F(l)co a * - v] 2 (k- 2 £ - l) + J( 0 } 


(124) 


where 


J(0 = 4m(oi a * - u^*) 2 J ^ F(S) dC - 4111(0^* - oy^ 2 . 


'“a* + 


F(lK* - ^ k _ 25 F(^) dS + 2m(u a * - a**) 2 f * k _2? 


d? + 


|l ( >r +1 Iq jW " “b *)* 5 ’ [?V + F(1) V - ‘Vjk 


-&■+ 


(“a* " t%*)k'^F(S)j^ [I(S)] n dS 


(125) 


Upon neglecting all terms containing A(£) , the first approximation 
to J(S) is 


J-(5.) = X J(?) = 0 (126) 
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while, neglecting terms containing the second and higher powers of A(£), 
the second approximation becomes 


%u> - (<%* - ^*)\w ♦ ( kV ^ - 

2 (<V* - a>b*)^D a * - ^-jk 2 " 2 ? 2 F(U 


(127) 


where 


•MU = ^(6) + 2 F(U d£ 

1 J 0 


5 k- 2 ^ 2 A(0 d£ 


Equations (124) for the pressure distribution, (101) for the veloc- 
ity distribution, and ( 103 ) for the temperature distribution in terms 
of the independent variable (; represent the solution to the problem, 
valid to terms in the square of the Knudsen number, provided the assump- 
tions of constant specific heat, of constant Prandtl number and coeffi- 
cient of viscosity, as well as of heat conductivity proportional to the 
first power of the absolute temperature, are valid. 

The solution (equation (124)) for the pressure is more simple than 
Schamberg' s solution in the sense that the present assumption u* = T* 
eliminates the explicit dependence of the integral in equation (ll4) 
upon T*. Figure 2 shows the dependence of 0 \}?( l) upon k as given 

by equation (123) for the values of the physical parameters selected for 
air in the section "Case of Air" and discussed there. 

The explicit determination of the expansions of the coefficient q 
and of the variables |, A(£), F(£), G(.£), and so forth in powers Of 
Knudsen number is carried out in appendix D. 


Jo(£) = 4kSn 2 


(128) 
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FRICTION COEFFICIENT, SLIP VELOCITY, HEAT TRANSFER, 
_ TEMPERATURE JUMP,AND PRESSURE RATIO 
Friction Coefficient 


It is convenient to define the friction coefficient Cf as 


C f - 


T r0 

-p U 2 
2 » u 


(129) 


The expression for is obtained from equation (50) which, on using 

equations (8l), becomes 



( 130 ) 


From equations (89) and (90) 


r*u* 


d _ q d 
dr* m d£ 


Therefore 


t i e u dco* 

T r0 - " 3T df 


where 




e 




va 


(131) 


(132) 


Substituting equation (131) into equation (129) and using equations (79), 
(80), and (8l), 


Re c f = - Hi - V 22£ 
e 1 km 


( 133 ) 
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where 


Re e = t] Re' 


va 


(134) 


From equations (C21), (C27) , and (C29) the value at the wall of the 
inner cylinder becomes 


Re e c fa 


_ 4[l + (jOfe* - io>b*)Kn e + 


k{l + k - (f< 1)/(1 " k 2} 

Similarly, at the wall of the outer cylinder 


(135) 


Re e c fb “ 


4k 

1 + (l m a* - !<%*) 

I^e + ( 2 V ' 

' 2Pb*) Kn e 2 

1 + k - 

£(!)/(! - k| 



( 136 ) 


Since there must be equal torque on both the inner and the outer cylin- 
ders, Cfb/Cf a = k 2 , as shown by equations (135) and ( 136 ). If 


C f 


0 c f I 


1 + 


l c f 


Kn, 


2 C f 



(137) 


then 



(138) 


Re e 0 C fb 


4k 

- [pF(l)/(l - k)] 


1 + k 


(139) 
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From equations/ 
the inner wall is 


Slip Velocity 

y (-8-1-)-,— (C28-)-,. -and_(_C29-) >_the s lip veloc ity at 


" iV + 2 u a* m l“a* “*« * 2%* (lto) 


Similarly, at the outer wall it is 


If = l u b* ^e + 2 u b* ^e 2 = (lV/^e + ( l41 ) 


where ^o^*, ^oi^*, 2 a) a * > ^b* are Siven by equations (C30) 

and (C31). 

Heat Transfer 

It is of interest to calculate the heat transfer between the rare- 
fied gas and the cylinders. If the dimensionless heat transfer q* is 
defined by 


3* = 


hq r 

^e^wa 


(142) 


where 


X = TlX 
e 1 wa 


(143) 


then from equations (53) and (81), 


3* = 


(1 - k)X* dT* 


kq dr* 

For Pr/c = Constant, X* = n*; so using equations ( 89 ) and (90), 


(144) 


3* = 


1 - k dT* 


kmr* d£ 


(145) 
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From equation (l45) and using equations (C22), (C27)> and (C29), the 
dimensionless heat transfer at the inner cylinder is 


1 - k 
km 


1 - W * 


(iV* - T vb * iV)Kn e * ( 2 T a * - T sb * 2 V)V 


o 1 - k 2 - 2m - G'(0) r / \ 

(7 - l)Pr M ^ 2 — —T— 1 + - 1 ^*) Kn e + 

[I-k2-F(l3 d L 


(l“a* - l^*) 2 ^e 2 + 2 (^a* ~ ZPb*) 1 ^^] 


(146) 


Similarly, that at the outer cylinder is 


q b * - 4i[l - T ub * + (jT a * - T wb * !T b *)Kn e * ( 2 T a * - T wb * 2 T b *)Kn e 2 * 


<7 ‘ ^ 1 ^ * 2 ( 1<% ‘ ' + 


If 


(l<%* - l“>b*) 2 ^e 2 + 2 (^a* ~ 2 a> b*) iai e 


q* = Q q*(l + x q* Kn e + 2 q* Kn e 2 ) 


(147) 


(148) 


then 


„ 1- k 

q * = 

0 a km 


Jl - T wb * + (7 - l)Pr wa % s 


1 - k - 2m - n G' (0) 


[7- k 2 - 0 f(i7] 


(149) 
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0 M b 


*“ = 



T^ b *-+- (-r - l)Pr wa My- a 2 



2k 2 m - qG * (l) 



k (15P) 


Temperature Jump 

From equation (C27), using- equations (C29) and (8l), the temperature 
jump at the inner cylinder is 


T - T 
■ l a ■'•wa 


■-wa 


= iV Kn e + 2 T a * Kn e £ 


(15D 


Similarly, the ,temperature jump at the outer cylinder is 


■^7"^ - lV + 2 V (152) 

wa 

The expressions for ]T a *> -^T^*, 2^a*’ an< * 2^b* are S ive n by equa- 

tions (C30) and (C31)» 


Pressure Ratio 

The ratio 'of the hydrostatic pressure at the wall of the inner 
cylinder to that at the outer cylinder is obtained from equation (11*0, 
. upon using equations (121) and (8l), as follows: 


Pb 

Pa 



exp 


^wa 2 ^ f 1 u*(cu *) 2 d£ 

^ J 0 T *(k2) 5+^(0 


(153) 
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For the case p x T, equations (124) and (8l) give equation (153) the 
form 


^ = [l + Kn e 2 0 *(lj] 


exp 




2rj|l - k 2 - F(l3‘ 


{&• 


- 4k 2 m) + 


2(1 - k 2 - 2m)F(l) + (k -2 - l)F 2 (iy] (^ a *) 2 + 

1^2(1 - k 2 ) + 2(1 + k 2 )m + (l - k -2 + 2m) F( l7] 20^*0^* + 
(k -2 - k 2 - 4m) 2 + J(l)} 


(154) 


In general J(l) is a small number and can be replaced by qJ( 1), which 
is given by equation (d4i). From equations (154) and (C29) 


o\ Pa / 




exp 


|(l - k^ - 4k 2 m) + 


2 olf - k 2 - o F < 1 0 : 

2(1 - k 2 - 2m) 0 F(1) + (k- 2 - 1) Q F 2 (1) + |j( lTj (155) 
For the case of equal wall temperature, T w v* = 1 and, from equa- 


tions (D23) and (C29), Q F(l) = 0* Therefore 


vm = 

. oV p J 


^W8 


exp 


2 Ml.- k 2 ) 2 


[(l - k 4 - 4k 2 m). + 2 J(1)J 


( 156 ) 


Or in terms of e 


:©■ 


™wi 


exp 


6 0 ^ 


( 1 + i 6+ i_ £ 2 + i t 3 a ( 4 + ...). 

V 2 io 5 7 / 


(7 - l)Pr 


120 


Wwa e / 7 lk 2 ? 2 \ 

r— ( 1 + I e+ 5T £ + • • •) 


0^ 


(156a) 
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From equations (C28) and (154), using equations ( 8 l), (C29), and 
(C27) and neglecting the first-order corrections to qF(i) and 
-it— follows— that 


lPb* = 


^ 1 - k£ 




1 - k 


^ - 4k 2 m) + 2(1 - k 2 - 


2 m) qF( l) + 


(k-2 - 1) Q F 2 (l)]( 1 a3 a * - i ^ + 1^2(1 - k 2 ) + 2(1 + k 2 ) 

(k-2 - 1 - 2m) oPdTJ - \ lTl * 2 J(l)j 


For the case of equal wall temperature Q F(l) = 0 and 


^a . 2 

0 n( 1 - k 2 ) 2 


- (l - k^ - 4k 2 m)^ 


lV - 2 l 1 **) + 


20 L + k 2 + (1 + k 2 )^'-^* - 2 j(i)J 


( 1 + i 6 + J.c 2 + I € 3 + l e 4 + 

3 o q |V 2 10 . 5 7 


• •)^l <u a* - | 1 ^*) 


(l + e + 3_ e 2 + it,3 + 5 ^ 


• lV 


€ m*(7 - l)Pr V /I 1 

Bo I \3 + 2 e 

0 ^ 


143 2 


( 158 a) 



NACA TN 2895 


h 8 


The use of effective values p e and X, e and hence of an effective 
Reynolds number Re e , similar to the effective Knudsen number Kh e 

defined by equation (109), is a natural consequence of the form of the 
expressions for the skin^friction coefficient, the velocity of slip, the 
heat transfer, and the temperature jump. The use of these effective 
values gives expressions closest in form to those of Schamberg and agrees 
with his when the functions F and G' are zero. Since q as deter- 
mined in appendix D will, for equal wall temperatures, increase with the 
Mach number and "curvature" l/k as shown in figure 3, so also 

will p e and X. e . 

The limiting values of all expressions when k = 1 or e 0 
agree exactly with the plane Couette flow solution of reference 7. 


CASE OF AIR 


Assuming that the gas is air and that the cylinders are made of 
metal such as aluminum or brass, the required physical constants are 
given as follows: . 

For the ratio of the specific heats take 

7 = 1.^00 (159) 


a value which is reasonably accurate and has the advantage of computa- 
tional simplicity. It also agrees with the value given by the kinetic 
theory for a diatomic gas whose molecules have five degrees 'of freedom. 

The kinetic theory gives, for the Prandtl number, Pr = ^-7/(97 - 5) = 
0.737 at 0° C for a diatomic gas with 7 = 1.400 (reference 4, p. 182). 
For air, various values have been used in the study of the laminar bound- 
ary layer in compressible flow by different investigators. The value 
Pr = 0.725 was used by Crocco and Conforto (references 15 and l6) in 
1941, Pr = 0.733* "by Brainerd and Emmons (references 17 and l8) in 1942, 
Pr = 0.750, by Schamberg (reference 3) in 1947; and Pr = 0.715; By Cope 
and Hartree (reference l4) in 1948. The value Pr = 1 was also used by 
several writers (references T9 to 2l) for mathematical simplicity. In 
the calculations to follow the value used will, be 


Pr = 0.715 (160) 

as suggested by Cope and Hartree in 1948, based on the latest data for 
air given by Kaye and Laby (reference 22), i94l. 
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Maxwell’s reflection coefficient a is given for air on machined 
brass by R. A. Millikan (see reference 4, p. 299) as 

a = 1.00 (l6l) 


The accommodation coefficient a has, according to M. Wiedmann 
(see reference 1, p. 658) the average value of 


a = 0.900 (162) 

for air on metal. This value is also used by Schamberg (reference 3). 

With these physical constants, equations (64) give the following 
constants for slip flow: 

a-L = 1.253 
c ^ = 2.498 

ej^ = 0.251 

eg = 1.911 

e Q = -5.491 

It is noted that, since air is not a monatomic gas, c-^' is used instead 
of c^. 

For the case that n ~ T, 

3 = 1 (164) 


(163) 


On using equations (C29), (159), and (160), equation (Dl6) gives 


0^ ’ K 1 + T „b*) 


. °*1^3 /l + k 2 l\ M 2 
1 - k 2 \l - k 2 m j wa 


(165) 
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For the case of equal wall temperatures , 


Twb* = 1 

6F(1) = 0 


(166) 


Similarly, using equations (C29), (159)> (l6o), (l64) to (l 66), and (D28), 
equation (D9) gives 




0.02045 

1 - 8k 2 + k^ 

9(1 + k 2 ) , 3 

(1 - k 2 ) 2 

2(1 - k 2 ) 2 

4(1. - k 2 )m m2 


(167) 


The variations 


of both -4ri and 

0 




for various values of the 


ratio k of the radii of the cylinders at various Mach numbers and 
equal wall temperatures with the ratio of the specific heats 7 = 1.400, 
Prandtl number Pr = 0.715, and the viscosity index 3=1 are given in 
figure 3. 


It is seen that for k > 0.40 and M < 4, ^ and have prac- 

tically the same value, and therefore equation (l65) can be used for 

^ O * , 

calculating q t) instead of equation (167). 


At smaller values of k and higher Mach number ^ ^ ; that is, 

the correction tends to decrease the value of i). 

From equation (156), on using equations (D4l), (159), (160), (l64), 
and ( 166) , 


log 


1 ^ 0-7000 A . k 4 . hA ) _ 

n h ^ / , ,9\2\ 


e ° ib 0 h (1 - k 2 ) 2 


Mwa 0.2002 
0 1 2 (1 - k 2 ) 2 !^ 


id + * 2 ) * -2^ - 

2 1 - k 2 m 


(168) 
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The- variations of ^p^* with the ratio k at various Mach num- 

bers and at equal wall temperatures with 7 = 1.400, Pr = 0.715, and 

3 = 1 are given in figure 4. It is seen that the pressure^ ratio ^p"”* 

Ob 

increases rapidly with the curvature l/k, especially at high Mach 
numbers . 

On using equations (C30), ( 166 ), and (l63), equation (l4o) gives 


l u a* 


'1.995 

k(l + k) 


(169) 


Similarly, equation (l4l) gives 


iV " 


l.995k 


(! + k) 0 p b * 


( 170 ) 


The variations of 1 u a * and 1 u b * with the ratio k are given in 

figures 5 and 6 , respectively. It is noted that u * is independent of 

the Mach number M , while u* depends on M through _p*. . As 

k decreases or the curvature increases, it is seen that u a * becomes 

more negative or the slip velocity at the inner cylinder ‘ increases in 
magnitude while the slip velocity at the outer cylinder decreases. 

Equations (C29), (166), (159), (l6o), (l63), (D31), and (C30) give 


l T a* = 


0.569 / 2 

• * 

• 


k(l + k)\i , k 2 


0.163 

1 + k 2 k 2 1 


k(l + k ) ( 1 - k 2 ) 

4(1 - k 2 )m ( X _ k 2) 2 2 m 2 

1 

o 11 


(17D 
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and 


rr> * 

Id 


0 . 569/1 2 k 2 \%s 


1 + k\m 


k 2 / oV 


0.163 

1 + k 2 k 2 

- -f- 

1 

v 1 * 

(1 + k ) ( 1 - 

k 2 ) 

4(1 - k 2 )m (1 - k 2 ) 2 

2mf_ 

■Hi p * 





0 l 0 P b 

The variations of 

l T a 

* and ,T. * with the 
1 1 b 

ratio 

k for 


(172) 


Mach numbers and equal wall temperatures with 7 = 1.400, Pr = 0.715, 
and p = 1 are given in figures 7 and 8 , respectively. It is seen that 
as k decreases pT a * becomes much more important than pTp*, the 

latter going to zero with k for all values of M. 


In figure 8 the anomalous behavior of the curves for small and large 
values of k is due to the large exponential values of qP-^* for 

increasing Mach number at all values of k except in the neighborhood 
of k = 1. In fact, for k = 1, pT b * = 0.2845M wa 2 , independent of Q p b * 

This is seen upon expanding equation (172) in powers of g = 1 - k 2 and 
letting g — > 0. 

From equation (DIO), the variations of ,q* with the ratio k for 

various Mach numbers and equal wall temperatures are found and plotted 
in figure 9. Since .^q* > 0, it follows from equations (C28) that q 

increases with the Khudsen number. 


Making use of equations (l)4l), (159), (l6o), and (l66), equa- 
tion (158) gives 


-n * 

l p b 


1 - k^ - It-k^m ^wa 7 \ 

0.700 — (2 pV - pH*) 


(1 - k 2 ) 2 0^ 


-1 + k g; + (1 + k 2 )m %a 2 
(1 - k 2 ) 2 o 11 

0.2002 ^r-7rli(l + k 2 ) + 


2.800 


l“b' 


* + 


(1 - k 2 ) 2 


2k 2 m 1 - k 2 

Mwa^ 

1 - k 2 m 

Qi 2 


(173) 
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The variations of n p * with the ratio k for various Mach num- 

1 d 

hers are given in figure 10. Since < 0, it follows that the 

pressure ratio P^*~ decreases as the raref action- of- the gas, -that -is,-- 
the Knudsen number, increases. 

From equations (135) to (139) and (l46) to (150), the variations 
of Re e 0 C fa and Re e 0 C fb , pCf* and jC^, 0 q a *, and 0 q b * with 

the ratio k are found and plotted in figures li, 12, 13 , and l4, 
respectively. 

* * 

The zero-order values of both the skin friction and the heat trans- 
fer are observed to increase in absolute magnitude at the inner and 
decrease at the outer cylinder with increasing curvature. 

Upon using equations (159)* (l6o), and (l66), equation (123) yields 


0^(1) - 


1-51 *vaT 
k 2 ( 1 + k) 2 


(i - k^) + [T- o(v)"5E 0 ‘ 136 ' _ 1,66015:2 


0.l69k^ + 0.06o('l - k 2 )(2 



(174) 


The variations of \|r(l) with the ratio k for various Mach num- 
bers and, equal wall temperatures are given in figure 2. It is seen that 
Qip(l) > 0, for sufficiently small k, and it increases rapidly as k 

decreases. It follows that the effect of the Burnett terms in the dif- 
ferential equation tends to increase the pressure ratio pjj*. 

The coefficient of the last term in equation (123) is plotted, 
against k in figure 15.” For the case of equal wall temperature, this 
term becomes zero. 

Consider in more detail the special case 


k = 0.5 


(175) 
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and 


p 

11 

ro 

From figure 3 

o n 

From figure k 

lp * 
Ob 

From figure 5 

l u a* 

hence 

l Q a* 

From figure 6 

l u b* 

hence 

1 <V 

From figure 7 

. l^-a* 

From figure 8 

l T b* 

From figure 9 

l 1 ** 

From figure 10 

. lPb* 

From figure 11 

1 

0 C fa Re e 


0 C fb Re e 

From figure 12 

l C fa 

and 

l C fb 

From figure 13. 

n * 
0 y a 

From figure ik 

0 q b* 

From figure 2 

0 ?(D 


1.17 
2.77 
-2.67 
-2.67 
0.24 
0.12 
3.73 
0.42 
0.96 
- 6. 30 

5.32 

1.33 
-2.79 
-2.79 
- 1.86 
0.59 
4.56 


> 
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( 176 ) 


(177) 
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From equations (C33), (C3b), and ( 113 ) 
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Equation (C32) gives then 

o\ - 5A6 
o \ - -°- 16 

’ o\ - 10 -° 
o^b '= -o.i4 

and equation (C3l) gives 

^ a * = 1T - 8 
2 0i b * = -0.08 

2 T a* = - 2 ^° 

2 T b * = 0.38 


(179) 


(180) 


From equations (D9)' to (Dll), using equations (177) and (180), 

q = 1.17(l + 0.96Kn e - 5.0Kn e 2 ) (l8l) 


Neglecting the rarefaction correction to \ and F(l), equations (135) 
and (136) give the friction coefficients 9 


Re e C fa = 5 v 32(l - 2.79Kh e + 17.9Kn e 2 ) 
Re e C fb = l«33(l - 2.79Kh e + 17.9Kn e 2 ) 



(182) 
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Similarly, equations (l4o) and (l4l) give the slip velocities 


- a Y~ U = -2.67Kn e + 17.8Kn ^ 

u e 


u h j O 

— = 0.24Kn - 0.l6Kn 2 

U e “ 


Equations (l46) and (147) give the dimensionless heat transfer 
q a * = 1.86(l + 10.5Kh e - 83.7Kn e 2 ) 

>■ 

= 0.59(1 - l-3Kh e + 8.2Kn e 2 ) 

Equations (151) and (152) give the temperature jumps 


— — = 3.73Kn - 27.0Kn 2 

m c e 

• - L wa 


— — = 0.42Kn e + 0.38Kn e 2 


Equations (C28) and (154) give the pressure ratio 


Pb 

Pa 


= 2.77(l - 6.3Kn e + . . .) 


The variations of q, Re e C f& and Re e C fb , (u a - u)/u 
u b /U, q a * and q„*, and (t, - T va )/ Tua and (t„ - T w „)/T va 
the effective Knudsen number Kn e for Pr = 0.715, 7 = i.4oo, 


(183) 


(184) 


(185) 


(186) 

and 
with 
3 = 1, 
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T wl3 * = 1, a = 0 . 900 , a = 1, k = 0.5, and = 2 are given in fig- 

ures 16 , 17, 18 , 19, and 20, respectively. Similar results could be 
calculated and plotted for any other choice of the parameters but the 
enormous amount of work required to do so for any range of parameters 
does not appear to be justified at present. This is especially true 
in that no experimental data exist for a check. 

The effect of the Burnett (second-order) terms or the Kn g ^ terms 

in equations (l8l) to ( 185 ) is seen to counteract the effect of the 
first-order slip terms or the terms in Kn e in all cases except the 
temperature jump at the outer wall. The latter would show the same 
behavior if and F(l) were not neglected. 


SUMMARY OF RESULTS 


A study was made to determine the effects of variable viscosity 
and thermal conductivity on the high-speed slip flow between concentric' 
cylinders. The results are summarized as follows: 

1. Satisfactory estimates of the effect of variable viscosity and 
thermal conductivity upon the velocity and temperature distributions 
were obtained from Schamberg’ s solution for constant values of these 
coefficients by basing the friction coefficient and the coefficient of 
heat transfer on the effective coefficients p e = T)|i wa and \ e = T|X^ a . 

These effective values of the coefficients p e and X e , in the case 

of equal wall temperatures, increased with the Mach number M and the 
"curvature" l/k. 

2. Only the expression for the pressure ratio p^p a was signifi- 
cantly different in this case from Schamberg’ s solution with constant 

H because of the use of effective values of Reynolds number and Knudsen 
number . 

3. The effect of the Burnett terms in the differential equation 
was more pronounced upon the pressure ratio p^P a . This effect 

increased with the effective Knudsen number Kh e , the Mach number M™, 
the temperature difference (T wa - T w ^)^T wa , and the curvature l/k; it 
was measured by the' factor ijf of equation ( 123 ) as used in equa- 

tion (12k) and plotted in figures 2 and 15. Effect of the Burnett 
terms could best be seen by an experimental determination of the pres- 
sure ratio p, /p . 

b / a i 
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4 . The effect of the Burnett terms both in the differential equa- 
tion and in the boundary conditions tended in all cases considered to 

counteract^bhe_firstj=.order_effect_of— s-l-i-p— ve-looi-t-y— and— temperatiire — 

jump on the boundary. 

5 . For equal wall temperatures the effect of the Burnett terms 
was to increase the pressure. ratio and to increase the skin friction 
on and the heat transfer to both cylinders. 

6. The curvature effect on the behavior of the friction coefficient, 
the slip velocity, and the temperature jump as k decreased was such 
that they all increased at the inner cylinder and decreased at the outer 
cylinder. When k approaches unity, the values of all quantities reduce 
to those for plane Couette flow. 


University of Washington 

Seattle, Wash., December 20, 1951 
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APPENDIX A 
SYMBOIS 


Adiabatic speed of sound, equation (34) 

- determinant of metric tensor, appendix B 
radius of inner cylinder 



A 



A i = 


F i 


A iJ 


numerical constant in boundary conditions, equations ( 64 ) 
fundamental metric tensor, appendix B 
constant of integration, equation (87) 

. . Ag numerical constants, equations (120) and (122) 

1 dp 
P Sx 1 

tensor 




nondivergent symmetrical tensor associated with A^, 
equation (21) 

radius of outer cylinder 

numerical constants in boundary conditions, equations ( 64 ) 
constant of integration, equation (87) 

l 

mean molecular speed, equation (38) 

numerical constants in boundary conditions, equations ( 64 ) 


C 


g C 0' g C l 


constant of integration, equation (9*0 
constants of integration 


Cf 


friction coefficient, equation (129) 
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D 

D/Dt 


6 ij 


E 

f 

F 


F i 


g 

G 

•h 

H 

I 

J 

k = a/b 
K lf Kg, . 


Kn 

Kn e 

Z 

L 


specific heat at constant pressure 

specific heat at constant volume 

constant of integration, equation (94) 

convective time derivative of hydrodynamics, 
equation (4) 

covariant derivative operator, appendix B 

• e po numerical constants in boundary conditions, 
equations (64) 

symmetric rate of deformation tensor, appendix B 

internal energy per unit mass of gas, equation (3) 
function, equation (98) 
function, equation (102) 

component of external force per unit mass, equation (2) 

function, equation (100) 

function, equation (104) 

gap between cylinders (b - a) 

enthalpy per unit mass of gas, equation (9) 

function, equation (D27) 

function, equation (125) 

■ numerical constants of stress tensor, equations (27) 
and (28) 

Knudsen number, equation (4l) (z/L) 

effective Knudsen number, equation (109) 

mean free path of gas, equation (39) 
characteristic length of flow 
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m = log e = log e equation (89) 


M 


Mach number 


N^, N^, numerical constants, equations (ll8) 


P ii 


Pr 


r q i 


(n) 


q i 


Q 

r 

R 

Re 

S 

t 

T 

Tw 

u, v, v 

u. 


U 


x, y, z 


X. 

1 


static pressure of gas, equation (5) 
component of stress tensor, equation (2) 

Prandtl number, equation ( 35 ) ( c pP/^) 

component of heat-flux vector, equation (3) 


rth-order approximation to q^ 


nth-order correction to ^q^, equation (l4) 

heat received per unit mass of gas, equation (8) 
radial distance from center of concentric cylinders 
gas constant 
Reynolds number 

entropy per unit mass of gas, equation (9) 
time, equation (l) 

absolute temperature of gas, equation (9.) 
absolute temperature of wall 

components of macroscopic velocity in x, y, and 
z direction, respectively 

component of macroscopic velocity, equation (l) 

surface velocity of inner cylinder (ao^) 

Cartesian coordinates of physical space 
component of Cartesian coordinate of physical space 
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X, Z 
a 


3 

7 


kZ 


ij 


functions, equations (75) to (78) 
accommodation coefficient, equation (58) 
viscosity index, equation (5 6) 
ratio of specific heats of gas ^Cp^c y ^ 
Christoffel symbol of second kind, appendix B 

thermal conductivity index 
unit tensor, equation (5) 

function, equation (93) 



r l 



0 


u 

P 


a ■ 

T 



equation (Dl8) 

constant, equation (9l) 

. 0c- numerical constants in heat-flux vector, equa- 
tions (29) and (30) 

divergence of velocity, appendix B 

coefficient of thermal conductivity, equation (20) 
* 

coefficient of viscosity, equation (19)- 
density of gas, equation (l) 

Maxwell* s reflection coefficient, equation (57) 
shearing stress 

component of viscous tensor, equation (5) 
rth-order approximation to t . . 



0 

0 


nth-order correction to ±y equation (13) 

angular coordinate in cylindrical coordinates 
dissipation function, equation (24) 
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V 


°Va 

cu* 

u3 

S' 

£ 


function, equation (115) 
function, equation (121) 
angular velocity of inner cylinder 

dimensionless angular velocity, equation (95) 
numerical constant, equation (C5) 
variable, equation ( 89 ) 
variable, equation (90) 


Subscripts: ' 

0, 1, 2, . . . . in front of any symbol, order of approximation to 

boundary conditions 

a, (3, i, j; k, Z, m covariant vector and tensor indices, appendix B. . 
e effective value 

wa. conditions at wall temperature T wa 

| indicates covariant differentiation, appendix B 

Superscripts: 

0, 1, 2, . . . in front of any symbol, order of approximation to solu- 
tion of differential equation, equation (110) 

a, p, i, j, k, Z, m contravariant vector and tensor indices, 

appendix B 

,p exponent in viscosity-temperature relation 
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APPENDIX B 

TRANSFORMATION OF EQUATIONS - OF - MOTION TO POLAR “COORDINATES 


The equations of fluid mechanics have been written in Cartesian 
tensor form in this report as in references 1 and 3. However, the 
problem considered had to be set up in cylindrical coordinates and 
Cartesian coordinates were abandoned in the process. In order to take 
into account all possible systems of curvilinear coordinates, it is 
best to express the equations of fluid mechanics, including the Burnett 
terms in general tensor form for any space with a metric form of the 
type (reference 23) 


ds 2 = a-y dx 1 dxJ 


By the principle of covariance all that one has to do is to express- 
each term in the equations as the proper invariant (scalar, vector, or 
tensor), which reduces to the known form when the coordinates are 
Cartesian. The method and notation used are to be found in reference 28 
The distinction between superscripts and subscripts is necessary. 

There is no difficulty with the equation of continuity. It is 
tim^derivative 01 ^ ° f “ SCal “ define the comoving 


D^ 

Dt 




where 

as p 
tive. 


D k is now the covariant derivative operator. For any scalar such 

the covariant derivative is the same as the ordinary partial deriva- 
The divergence of the velocity is the scalar 


1 


0 = D u k = u k 

k |k 



( 
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In other words, take the covariant derivative of the contravariant 
velocity vector u k and contract the result; a = | a i j j ^ ekerm i _ 

nant of the metric tensor By definition the covariant derivative 

is 



6u k + p k 
3x Z l] 


U 


m 


• 

Where the terms r , are the Christoff el symbols of the second kind. 

Zm 

The continuity equation is then 


or 


Dp 

Dt 


+ p© = 0 




(Bl) 


for, when covariant differentiation becomes ordinary differentiation, 
these expressions reduce to the known equation, equation (l) (refer- 
ence 23 , p. 196). 

Similarly, the momentum equation is a vector equation and the energy 
equation is a scalar equation. All that is needed to give them general 
tensor form is to place the indices properly and to use covariant dif- 
ferentiation. They are obviously 



p F k + 


8x k 



(B2) 


DE 

p — + p0 
Dt 


cx3 

+ T U 


a|p 


+ q 


a 


a 


= 0 


(B3) 


which replace equations (6) and (7). 
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Expansions (13) and (l4) remain the same as long as it is borne in 
-mind that e ach j ( n ) is a covariant tensor and n ) is a covariant 

vector. Equations (19) and (20) for t lie s econd~order~ appr ox-i-mat ions 

remain the same: 


i T u - -j*i 


8xJ 


iQi - -X- — r 
Sx 1 

8ui 

provided ^ — is redefined as the covariant tensor 

3lt j 


^Ui _ 1/ \ 1 

a x . “ e± i 2V 1 u + u ji v ~ 3 0a 

J 


3 ^ 


One must use a.j_j here as it is the fundamental tensor of the space; 
a 1 J is the fundamental contravariant tensor, and 5^ is the funda- 
mental mixed tensor. 

By continuing this process the expressions for the third-order 
Burnett terms become 


T ij (2) = K 1 7 0e ij + K 2 ! 7||( A i|j + A J|i) - f A a |a a ij - |( u “|i u j|m + « m |j u i|m) 
I uk !m ““Ik a ij - ( e i m u m| j + e j m u mli) + § ^ u klm a J + 

K3 ^uj + K ^|l( p ii t u +p u T ii) - 

3 ^ P l* "I- + *3 ^( T U T l0 * 5 akm *|k V a ij) + 

|( e i a e aj + e j a e ai) ' 3 e “ P e aP a ij] 


k 6 u_ 


q . ( 2) = e i 0T li - e 2 J^er) + 2u* (i t,^| + e 3 g P)k 


e k . 
e i + 


~ ek i|k + e 5 T |k ek i 
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where the covariant vector 
and 



1 6p 

P cbc 1 ' 


m 8T 8p 

T|1 “ a?’ P|i ' 


= 5 s ! r k 6r 

^ Sx 1 8x k 


The Cartesian symbol A^ * has been replaced by the covariant expression 


1, \ 1 
L . = -(A. . + A.. ) - - 
iJ 2\iJ Ji/ 3 


A . a. . 
3 k iJ 


In order to obtain equations (5^-) to- (56), it is necessary to intro- 
duce plane polar coordinates r = and 0 = x^. The metric form is 


ds^ = dr^ + r^ d0^ 


11 2 1 

so a-Q = a = 1, a22 = r = — and all other values of a^ are 


p 

zero while a = r . The only nonvanishing Christoff el symbols are 


r 22 = " r 


_2 n 2 1 

r 12 r 21 r 


Denote the physical components of the velocity vector by u r and Uj, 

o Up ^ 

* = Similarly, the physical components 


where Uj. = u^_ = u^, u 0 = ru 


of the vector are q r 

ence 23, pp. 1^5 to l46). 


q0 and, of the stress tensor (refer- 


and 
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rr 


= T 


11 


.11 



12 _ 



Then 


00 


-2 T 22 


2 22 

rr 



au r | u r ( 1 c> u 0 
Sr r „ r S0 


and the equation of continuity becomes 



Written out, equations (B2) and (B3) become equations ( 31 ) 
In covariant form D k = a ka D a , 


P 


Du k 

Dt 


pF k + 


Ik 



+ D a T 

ka 


= 0 


* DE 0 aB 

P frr + P© - 2)ae p u 
Dt P|a 


- c >v) + 


.(2)ap 


u 


.( 2 ) 


a 


a|0 


l a 


= C 


Du k 

The components of become, 

Dt 

nonvanishing Christoff el symbols. 


(B4) 

(B5) 

and .( 32 ) 

(b6) 

(B7) 


using covariant derivatives and the 
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DU]_ _ 8^ k _ 8uj. 8u r u£ 8u r u0 2 

Dt 8t U U l|k 8t r 8r r ^0 r 


Du„ 8u 


8u„ 


Dt 8t 


* sr + uk "sin * r sr + "r + y 


buy 

80" 


Let F x = F r , F 2 = rF^, D-jp = gL, and D 2 p 
ponents of e^j are 

8u r 1 . 

611 ' ST " 3 


The covariant com- 

80 


®12 e 21 


l/8ur 8u0 \ 

= + r 87 ' u t) 




while the contra.variant ones are 


11 _.^ u r 1 

6 " Sr " 3 


.12 


,22 


21 _ l/ 1 ^ u r 

1 Su0 

• 2^-2 ^0 

r 8r 


l 1 1 

' ?fer “r; 

1 ' 3? 


Since 


du T 




u 


111 ' 8r U l|2 80 " U 0 


8u0 


U 2U = r &T 


/8u0 \ 

U 2I 2 = T \W + Ur J 
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it follows that in equation (B7) the dissipation term is 



V) 2 + (1 

r / \r 90 r ) 


1 

3 


0 


2 


Consequently, the equations of motion’ and energy in two-dimensional 
polar coordinates are 


duy 3u r U0 8uj. 

,3t ^ 3r r 80 






18 fl 8uj. + Su0 
r 8^ \r 80 8r 



2pp u r 
r \8r 


1 ^0 
r 80 



■+ D 


a 



= 0 


(B8) 



P F 0 + 


1 8p 

r 80 


iQi ^ + 5a _ ^ 

8r \ r 80 8r r / 


.a i Ji 52 + 2 e . i 

r 80 \ r 80 r 3 / 


2p/l ^r + ^0 
r \ r 80 8r 



+ D 


a 



= 0 


(B9) 


/8e 8e u 0 3 e\ 

( — + u- — + — — V + 

\3t 8r r 30/ 


P0 



l/l ^r + ^0 _ ^0\ 2 + 

2\r 30 3r r / 



I A( rX ' 
r 3rl rX 37 J 


r^ 30 \ 30/ 


r ( 2)kZ 


u 


k | 2 


J2)k 
q Ik 


= 0 


(BIO) 


For 0=0 and p = Constant these reduce to the well-known equations 
for an incompressible fluid (references 12 and 13). It still remains 
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to find the components of and t^ 2 ) a.nd substitute into these 

equations. For the former the result is immediate: 


q (2) - e >£e^ + e v £ 

r 1 pT o r ^ pT 


!*“*> 


Su r - dT 2 
Sr Sr + 


^ ST 

r Sr S0 


-1 





(2) 


f) ti! n i ® + fl P 

1 pip “ ?\(A 2 — 


1 St 

r 80 


■2 


I i A( 0T ) + 

3 r 80 



r / Sr 



+ 0, ^ 
“* PP 



^ u 0 U0\8p 

Sr r /Sr 


.L 0 + (i ^ + iiU ip 
3r S0 \r S0 r / r 80 


+ 01. !=- 


1 j^Yl + ^0 V| 

2 Sr \r 80 Sr r / 




+ 0 C 


pT 


l/l ^r 

2\r S0 


+ 


^ u 0 

Sr 
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A longer, more tedious calculation leads to the stress components 


T 


( 2 ) 

rr 




f i ^ + M 2 + ifi H _ y \ 2 

\r d0 r / 3^r d0 " r / 


2 / au r\ 2 _ 2 ^0/l Su r u 0\ 2 ( dn r 

\Sr / 3 dr \r d0 r ) + 3 ® \dr 





12 

-9 





+ K 0 


ti£ nfe s2 T 

P V dr 2 


_i_ dr 

3r dr 


1 + ki 2 / 2 ap &r 1 a P ^\ 

3r 2 d0 2 / ^ PPT\3 dr dr ^ p 2 5^ d0y + 



fdr \ 2 

^dr/ 
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Since for the symmetrical case considered in this report u r = 0 
and all partial derivatives with respect to 0 vanish, 0 = 0, 

A, = - i and Ap = 0, it follows that these expressions simplify 
■ p dr 

considerably, becoming simply 

q r < 2 > =0 

- 26 , id 'l ® + 1 , £ ivM _ ^ + 1 g k £ ±hi . y\ + 

. pT r dr 2 J pp dr\dr r J 2 ^ p dr\dr r/ 

e .i]a, d 

A pr \dr r J 2 5 prp dr \ir r / 
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Equations (B8), (B9), and (BIO) reduce for this symmetric case to 



l 0) ( 2 ) 

-)’ T 00 = T ff , 

case^ so these are pre- 

(B12) and (Bll) agree 
r terms occur in equa- 
tions (50) and (53) since r^) = 0 and (2) = The second _ 

( 2 ) ' 

order term for q^ does not occur because its derivative with 
respect to 0 is zero in the problem considered. 


But T = T (2 > T<l( - '(D / dU 0 U 

^r - = d p-» anci ^0 = Q0^ 2 ^ in this 

cisely equations (47) to (49) while equations 
with equations (51) and (52). No second-orde 
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APPENDIX C 

REDUCTION OF BOUNDARY CONDITIONS TO DIMENSIONLESS FORM 


Making use of equations (79) to ( 8 l), the boundary conditions, 
equations ( 71 ) to ( 74 ), can be written in dimensionless form as follows: 


V = 1 + a i 7 — o7i 


(t *' 


a ) 


d /u 
dr* \r 


- (‘ 

O) 4* — 


*>|2 
9 / v 


CJa. 


M' 


X„*d5^ (Cl) 


u b 


* = 0 - a. 


(V) V2 




T a *Pb* 


r* 


• d / u* 
dr*\r*y 


S + ^3 - g X b *E 2 

b (v) (v) 


(C2) 


V * 1 


a c. JtsLJSl) S + Z a *5 2 .■ 
(v) 1/2 V dr *4 (T a ‘f 


(C3) 


M V V /aT»\ E + h*) : 

T a *Pb* W*/ b (V ) 2 ^*) 1 


T * = T -u* - c n v — - ^ — 1^- ) d) + Z, *0)2 (C4) 

b wb 1 m /__ ^\2 /m * \2 


where 


(JO = 


= (1 - kb 1 / 2 

a Pa( T ya ) 1/2 k Bewa 


- °- 796 ( i r ± ) 


Kn wa ( c 5) 


using equation (4l). 
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The functions X a *, X b *, and Z^* are obtained from equa- 

tions (75) to (78) by. use of equations (80) and (8l) and are 
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• 7 8 


On using equations ( 89 ), (90), (91), and (C5), equations (Cl) 
to (C4) become , respectively, 


• u a* = 1 + 


0 . 796 a 1 (l - k) / dtD *\ 


kia^Ta *) 1 ^ 2 0 


(— ) + x a ^e 


(C7) 


V = 0 - 


0 . 796 a 1 (l - k)(T b *) 

M a *Pb* 


1/2 

/doo*\ 

tel 


Kn e + X b Kh e e 


(c8) 


T a ^ = 1 + 


0.796c 1 (1 - k) 


Vd^ /q 


(—) Kn e + Z a Kn e 2 


(C9) 


m * m * 

x b 1 wb 


Pb 


(V) 

' 1/2 I 

^dT*\ 

1 

* 

^ K 


Kn e + ^Kn^ (CIO) 


where 


^e = 


(Cll) 
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and 


0-6336(1 - k) 2 (n a *) 2 x a 


kV(T a *) S 


0.6336(1 - k) 2 
k 2 m 2 (T a *) 2 


: i"^-(^-i*)Sri^ 3 *Sr 5 ( ^.**>- 


"Ha' 

n 


'.(la* *21 + 1 Wl-_ 5 H 2 m< toA 
\8 d$ 8d( 6 7 ^ a d£ ) 


115=0 


_ 0.6336(1 - k) 2 (^ b *) 2 X b 

k2r > 2 ( T a*) 2 (Pb*) 2 


_ 0.6336(1 - k) 2 5 „„ d 2 co* 

- - 7T L* 

tan 2 ( T a *) 2 (p b *) 2 L 


(‘- 631 - i »)ff f - 5 T - If if 


P*) - 


dT * 


8i) at, 


-k= 1 


^ - 0 -6336(l - k) 2 (u a *f Z a * 

k 2 n 2 (V ) 2 


0.6336(1 - k)^ 
k?m2(T a ^) 2 


M ^ /dT*> 2 

e n 7M TQ T , „*( + eo ( 

r va a U I o\ 


14 d^ 2 


n 


ft '"-A*"* if * (f * *6>* if ; *4 * f "*)? 


5=0 


* 0.6336(1 - k) 2 (|i b *) 2 Z b * 

= k V(T a f (p,*) 2 


_ 0.6336(1 - k) 2 


k 2 “ 2 (v ) 2 (Pb*) 2 



k 2 e 8 f^f . d k^T* ^ 
»VdC ) 14 d ^2 


*1 


W ^ ^ egjkV 


5-i 


because of equations (C 6), (89), (90), (91), and (95). 
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>(C12) 
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From equations (102) and (98) 


r;t, 

fer 


(C13) 


F"U) = V” 1 [a( (f)] n_1 Qi A'(C) - 2m A(C[] (Cl4) 

n^I nl 


where the prime denotes the derivative with respect to Equation (93) 

gives 


A' ( 0 = — - 1 


(C15) 


Therefore 


F'(0) = 0 

F ' ( l) = 0 


(C16) 


F" ( 0) = 4m 2 - lj 

F" ( l) = 4k 2 m 2 ^- - 


and from equations (104) and (100) 


G* (0) = g'(0) - g(l) + g(0) 


(C17) 


= 2mf(0 - g(l) + g( 0) 


(C18) 


G"(£) = 2mf'(C) 


(C19) 


whence 


G" ( 0) = G" ( l) = 0 


(C20) 



NACA TN 2895 


81 


From equation (101), using equations (Cl6) and (C17), 

(-2m) 



-<»£*— mg* 



d^CD* 


^d 2 ay* 


1 - k 2 

- F(l) 

“a* - 

V 

1 - k 2 

- F(l) 

CD * 
SL 

- V 

1 - k 2 

- f( 1 ) 

“a* ■ 

' “b* 


( -2k 2 m) 


(**“ 2 H a */ri) 


(kk 2 m 2 \±^*l T\j 




\d£ 2 ) 1 1 - k 2 - f( 1 ) 

while from equation ( 103 ), using equation (C20), 


/dT*\ 

W h 


SSI) = - (t„* - T- 
d£ Jo { 


a ) 


*\ _ (y - l)Pr wa. M wa 2 (V - %*)‘ 


(l- k 2 -F(lT) : 


[l -k 2 - 2m -G' (oj 



= - (t & * - T b *) - - 7 ~ 1 - Pl ’ wa Mwa ' 2 ^ a * a>b *) 2 _ 

1 [I - k 2 - F( lj] 2 


k 2 - 2k 2 m - G’ ( lT] 


'd^A ( y - 1 ) Pr va\a 2 (V-V) J 


^£ 2 /0 


|l - k 2 - F( lT| 


(4m 2 ) 


/d 2 T*) 

| _ (’'- 1 > 1>r »a 

^K* - <%*Y (^V) 

W! 

1 

[l - k 2 - F(lj] 

1 2 


(C21) 


(C22) 
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Putting equations (C2l) and (C22) into equations (C7) "to (CIO) 
yields 


<v - «>b* 1 - k 


ao a * = 1 - 1.592a-. 


(V) l/2 k[l-k 2 -F(l[l 


Kn e + X a Kn e 2 (C23) 


= 0 + 1.592a-^ 


(“a* - k 2 (l - 


k) 


Kh e + kX b Kn a 2 (C2U) 


T a* p b* 1 - k 2 - F(l) 


T a * = 1 - 0.796 Cl 



T * 
a. 


V + 


, 0 1 - k 2 - 2m - G'(0) 1 0 

(7 - l)Pr wa M^a*- a**) 2 — 7 r + Z a Kn e 


[I - k 2 - F(lJ) : 


(C25) 


Tv* '= T wb * 


, (1 - ^(vO 1 / 2 

+ 0.796c, < 

-L wJ""P -X- "X - 


n£T a * P b * 


V - T b* + 


(7 - DPr wa Mwa 2 ^* - 


1 - k 2 - 2k 2 m - G' (1) 


[l - k 2 - F( lj] 


> Kn e + ^ Kn e 


(C26) 
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It is convenient to express o> a *, ' T a *, u u * 

Pb*, and tj in ascending powers of the effective Knudsen number Kn . 
Thus write e 


^a* = 


= (fa^ 1 * 1<V* K“e * 2 V ^e 2 ) 


<u b* = 0%* + i“b* Kh e + Kn e 2 


T& * = 0 T a^(l + lT a * Kn e + 2 T a * Kn e 2 ) 
T b* = 0 T b* (i + l T b* + 2 T b* Kn e 2 ) 


>■ (C27) 


u a* ~ O u a* ( x +. l u a* Kn e + 2 u a * Kn e 2 ) 


U *>* = O 1 ^* + iV ^e + 2 V 


P b* ~ QPb*^ + l^b* Kh e + 2 I>b* Kn e 2 ) 


> (C28) 


11 = O^i 1 + 1^* ^e + 2^* ^e 2 ) J 

where the subscript 0, 1, or 2 written in front of a symbol denotes 

the order of the approximation to the boundary values. 



84 NACA TN 2895 * 

Equations (C23) to (C26) give in successive steps the zero-order 
approx imat ion : 


* = 1 

0 a 


cA 


* = 0 


T * = 1 
0 a 


rp * = T * 
Oh wb 


(C29) 


the first-order correction terms: 


l^a* ~ “ 


1.592a 1 (l- k) 
k[l -k 2 - 0 F(lT] 


1%' 


1.592a 1 (l- klk 2 ^*) 1 / 2 
[I - k 2 - 0 F( l]j 0 Pb* ' 


, 1 - - 2m- qG t ( 0) 


0 • 796 c -1 ( 1 — k) j v, p ■ 

,T * <1-T wb +.(7-l)Pr W aMwa 2 - ; 

la km 1 [l - k 2 - 0 F( 1_)J 


HC30) 


l T b* = 


o. 796 c 1 (i- k) 


r 


m nP' 


Ob 


( T wP) 


,\!/2 


Jl-T w ^+(7 - 1 ) Pr wa Mwc 


3 1- k 2 - 2k 2 m- 0 G’ (1) 

[I - k 2 - 0 f( iT] 2 
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where qF(1), qG'(O), and qG’ ( l) are obtained from equations (102) 

and (Cl8) using the zero-order approximation, equation (C29), to deter- 
— mine— r \ — and and the second-order., correction terms _ 


2 a a* = 


1.592 ai (l-k) , 1 \ 

( -.co * - ox * - ^ _T *) + > 

k|I-k 2 - 0 F(iTl' a 21al 0 


1.5928^1 -k)k 2 (T ub * 1/2 

£“b* - — r - g Til (lV ' l“b* - l T a* + 

E- k2 -o F <^ oV V 


§ l T b* - lPb*) + * o x b 


T * = - 

2 a 


0.796c 1 (l-k) [/, , \ 

— to— jU 1V-21V w- w ivj 


( 7 - l)Pr 


1- k 2 - 2m- 0 G’ (0) 


[I-k2- 0 FUT] 


2 ,a>.»-2 1 V -5 iV) 


2 l la 


>(C31) 


V + n Z : 


0 a 


0.796c,(l-k) I • ' 

2V = TTiVwb iV + M 1 " 3T wb")lV - 

“ .Mv) L 


(* - T wb*) 


iV + (7- 


l)Pr M 


2 1 - k 2 - 2k 2 m - 0 G’ (1) 

(I-kS-oFdT) 2 


( 2 iV - 


2 iV* 5 iV-iV 



O^b 

^wb* 
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From equations (012) and (C29), with ti* = T*P, 
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APPENDIX D 

DETERMINATION OF APPROXIMATE EXPRESSIONS FOR DISTRIBUTIONS 
From equations (90) and (91) 



(Dl) 



(D2) 


The kinetic theory of gases (reference 4) gives 

H « T 13 I 


that is, 


y 

ii* = t*3 


(D3) 


Experiment shows that for most gases the value of (3 lies between 0.5 
and 1, and it varies slightly with temperature. Substituting equation (D3) 
into equations (Dl), (D2), and (93) yields 




NACA TN 2895 


89 


Q 


It is customary to assume (3 constant and for most gases its value is 
quite close to unity. While the calculations can be carried out for any 

constant value _of_ 0 _,_ it is . much simpler to take (3 = 1 as is done in 

the text from equation ( 124) on. " • 

Substituting equation (103) into equation (D5) and taking 0=1, 


•ft = 


_ (7 - l)Pr M 2 /o) * * ox*V 

V?-|(V-T b .) E 2 + _ ^ 1 


[l - k 2 -,F(lT] 


[c - |<i - 


* 2 >t 2 - ^ 


+ KO 


(D7) 


'5 


where l(£) = /' G( 0 d£. 

J 0 

Substituting equation (103) into equation (d 4), with 0=1, gives 


" ' f( T a* + V) * (7 - DPr ua -i 




%*) 2 IT 


[l - k 2 „- F(lJJ 2 l£ 


-U + k 2 ) - 


1 - k 2 
2m 


+ KD 


(D8) 


On neglecting the rarefaction correction to F(l) and l(l) and using 
equations (C27), (C28), and (C29), this yields as successive approxima- 
tions to rj 



+ 


(7 - l)Pr wa 

E - k2 - o f(i >] 2 



+ 


k 2 ) 



0 



(D9) 
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(7 - l)Pr Myg 


O' 


- §(lV + A*) * 2 (iV - A‘) E . k2 _ ^ 


|(1 + k 2 ) 


1 - k c 
2m - 


+ nKD 


(DIO) 


0^ 2^ 




T * + T * l + 
2 a 2 b 1 


') + (l“a* “ lV) + 


2 (2<%* - 2%*) 


(7 - l)Pr % £ 


- 1 [l - k 2 - q F( lYJ 


1(1 + k 2 ) _ l_il + 1(1) 

2 2m 0 


(Dll) 

Now A(0) = A(l) =0 and , for 0 < £ < 1, |A(£)| « 1. As a 

first approximation assume (superscripts denote order of approximation 
to solutions of the differential equations regardless of the boundary 
conditions) 


A(5) = X A(U = 0 


(D12) 


It follows from equations ( 98 ); (100), (102), and (104) that 


!f(S) = x g(5) = ^(S) = htt) = 0 (D 13 ) 


Substituting these into equations (101), (103), (DT), and (d 8) gives 



(Dlk) 
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V* = T„* - | 

(t * - T-l,*' 

L— a _h . i 

(7 - DPr My a 2 1 

K* - %*! 

ifl 

1 - (l-k 2 K - k 2 ^] 

OJ 

CVJ 

I 

1 — 1 
cu 

1 



(D15) 


1_ _ l/ m „ ,\ (7 - DPr M^ 2 (o3 a ^ - u^*) 2 


S = |(T a * + T b *) + 


2(1 - k 2 ) 2 


1 + - 


2 1 - k c 


m 


(di6) 


i, = 1( V ♦ V ) + ^ _ l)Pr ^( V . % .) 2 A * i 


€ + 


i£ e 2 + ± e 3 + 

60 ko ' ' * 


(D 17 ) 


where 


e = 1 - k £ 


and 


X| = f-K't; - 5(1,* - v)t 2 + 


(Dl8) 


(y - DPr Mwa 2 ^* - “>b *) 2 


i - |(1 - k 2 )C 2 


(1 - k 2)2 

For the second approximation, equation (d 6) gives 


1 - k 2 ^] 

2 m J 


(D19) 
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which, upon using equations (D15) and (Dl6), yields for 3-1 


2 a(0 =— 1 
2 Iq 


(v - V )^ 1 - 0 + 


(7 - l)Pr M wa 2 (c^ l * - (a^y 


(1 - k 2 ) 


Cd - 5) - 


1 ■- £■ + k 2 £ - k 2 ^ 

(l - k 2 )m _J 


(D20) 


where lq is given by equation (Dl6). It is noted that 2 A(0) - 2 A(l) - 0. 
Neglecting the second and higher powers of A(5), equation ( 98 ) yields 


2 f(£) = 4m 2 / k 2 ^ 2 A .(£) d£ 


which, upon using equation (D20), gives 


k 2 ^ , 

2 f(5) = f-4 


(T a * - T b *) + (1 - m)g + m£ 2 

' Pm 


(7 - l)Pr - q>t*) : 


(1 - k 2 ) 


-m£(l - 0 + 


1 + k 2 - k 2 ^ 
2(1 - k 2 ) 


Constant 


(D21) 
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It follows from equation (102) that 


2 F(0 = 


(V - V)[|d - rn)(l - k 2 0 + (1 - m)£k 2 ^ + 2 k 2 ^j + 


(7 - l)Pr % a 2 l> 


2 (m * - n\ 2 


^a* - “b*) : 


1 - k £ 


■*5(1 - Ok 


2£ 


(l - k 2 ^ )(k 2 - k 2 ^) 
2(1 - k) 2 


and 


(D22) 


T(l) • -l-(T a 

2 V- 


•)(■ 


* - Th* 111 + k 2 - 1 ' ^ 


m 


(D23) 


F(1) v) ( 1 + i , + i* 2 *w* 3 * - • •) 


(D2M) 


where e = 1 - k^ as before. 


Substituting equation (D2l) into equation (iOO) yields for 
Pr/cp = Constant 


2 g(C) = 2m / 2 f(0 a£ 


k£ 

V 


"(V - T b *) 1 


LJ5 + (2 . n ; h + . ^_ 1 . )Er - ^) 2 


2m 


2 


1 - k* 


1 - k 2 + 2k^m 
_2(1 - k 2 )m 


+ (l - m).£ + m£ c 


k 2 S 

Ml - k 2 ) 


+ g C 0 + gPl 5 


(D25) 
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From equations (D25) and (104) 


' G(U = ' V)j- 


3(1 - k 2 ) 


1 + k 2 k 


- lj (l - k 2 ^) - (2 - mKk 2 ^ - m^k 2 ^” 


(y - l)Pr My^K* - %*) 2 / 1 - k 2 + 1 + k 2 V 
1 - k 2 \ 2m 4 / 


— + ■ ■■ — 

1 - k 2 


m£ 2 k 2 ^ - - 1 - k 

4(1 - k 2 ) 


1 - k 2 ^) - (l - m)^k 2 ^ - 


(D26) 


The integrals 1(0 defined in equation (D7) for £ and l(l) in 
equation (d 8) for q can now be found from equation (D26). They are 


1(0 = / 2 G(0 dS 

J o 




f • iy - 


A(l . 1)(! . *25) . I 

2m \2 m J\ / 2 


1 - * i i 2 * 25 


(y - l)Pr Mw^V - o>b*) 2 J 1 - 4k 2 + _1_ 
^■T] ( 1 - k^) j _ 4(1 - k^) 2m 


fl + k 2 1 - k 2 

{ 8 Urn 




1 _ r 

2 m. 


+ | £ 2 k 2 ^ + ■ - 1 - k ^ 

2 16(1 - k 2 )m 


(D27) 
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KD = f 2 G(0 d £ 

J 0 


- - V)^! - * 2 > * * 

(r - Dpt . ^[± 1 8kg - t4 + _ j“| , 

U3(l - k 2 ) 2 16(1 - k 2 )m km 2 

J 


(D 28 ) 


KD = - 


120 ^ 


V ' MC 1 + e + • • •) 


(7 - l)Pr M^ 2 ^* - o^*) 2 (o + -i- e + . . . 


(D 29 ) 


The derivatives G'( 0 ) and G'(l) needed in the coefficients X T *, 
l T b*> 2 T a*' and 2 T b* in boundary conditions (C 27 ) are also found from 
equation (D 26 ) as follows: 


-G'U) 


— ^(T a * - T b *)[l + k 2 - 


2 _ 3 ^ 1 ~ *-l + (1 - m)k 2 ^ + 
2m 


2 m ( 1 - m) £k 2 ^ + 2 m 2 (; 2 k 2 ^| + 


(7 - l)Pr N^a 2 ^^* - c%*) 2 |"l' + k 2 1 - k 2 1 + k 2 


1 - k £ 


1 - k £ 


2 m 2 £( 1 - Uk 2 ^ 

1 - k 2 


(D30) 
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Thus 


2 G'(0) = 


(V - V) 


2 + k m - 


3(1 - k 2 ) 


2m 


(7 - l)Pr M^ 2 ^* - 0 ^*)' 


1 + k 2 k 2 m 1 

4(1 - k 2 ) (l - k 2 ) 2 


2 G'(0) = 


24 \ 


( T a‘ - V)(l> 15 ‘ + I * • • •) 


jo ( ' r - 1>Pr “ua 2 ^* - ''ll') 2 (i + f 6 + • • •)] 


and 


2g’(1) = — 1 

i T| 


(t 9 * - T b *) 1 + 2k 2 + k 2 m - — — ^ k ^j + 


(7 - l)Pr M wa 2 (m a * - a^*)' 


1 + k‘ 


k^m 


4(1 - k 2 ) (1 - k 2 ) 2 


2 G’ (1) = — — |7t * - T *^(l + 2- 
24 \ 11 a • b A 10 


, + f , 2 + 

4 


...) 


jtfr - Mwa 2 K* - V) 2 (i + | e * • - •)] 


(D31) 


(D32) 


(D33) 


(D34) 
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Using equation (D22) for 2 F(£) and equation (D20) for 2 A(£), 
integrals and J 2 in equations. ( 128 ) become 

•MS) = 2 F(S) + 2 F(C) d£ 

J 0 




(V - v)| 2 (i - »)S + 1(3 - |)(i - k 2 ^) 


(3 - m){;k 2 £ + 



(7 - PPr Mwa 2 K* - %*) 2 


1 - k £ 


2k 2 n 


1-1 


^•S^ t - k2£ ) + " jVt+(2 -- )Sk2 ' 


Jt (l) = - ■ 7 — 


( T a* ' T b*) 


1 + 1 k 2 _ ^ ^ + 

2 2 2m 




(7 - l)Pr M^ 2 (cu a * - CDh*) : 


1 + k 2 2k 2 m 

4 - 


2(1 - k 2 ) (1 - k 2 ) 2 m 


Jl(1) ' nt li v ' v) ( x ♦' 1 • * S ' 2 * S « 3 ♦•••) 


(7 - l)Pr - ®b*) 




the 


-i- 

.2 

(D35) 


(D36) 


(D37) 
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(2 + m)£k“ 2 ^ - m^k"^ 




(D38) 



(7 - l)Pr 



2 



+ 


2 k 2 m 

(1 - k 2 ) 2 



(D39) 



(7 - l)Pr - <v) 2 ±-(l + ± c + 



(d4o) 
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From equation (127).> when use is made of equations (C29)> 


0 J <?) = 0 j l ( ?> * 0 J 2<?> - 2k2 ' 2? o F(?) 


and so 

0 J(1) ■ 0 J ! (1) * 0 J 2 (1) - 2 O^ 1 ’ 


0^ 


( ^ ^wb*) 


9 + 

2 


2k^ - ir - 2m - 


(1 - k 2 )(7 - k 2 ) 


2 m 


)_ 


(7 - l)Pr % 2 


|d + 


k 2 ) 


2k 2 m 





(d4i) 


(d42) 
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Figure 1.- Coordinate systems. 
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Figure 2.- Burnett terms correction to pressure ratio 
diameter ratio k for various values of Mach number 

T wb * = 1$ a = 0.900$ 0 


Pr = 0.715$ 7 = l.kOOj 3=1$ 


p^* against 
Mwa* 

= 1 . 


gp ^3 
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Figure 3 .- Variations of ^tj and qT| with diameter ratio k for 
various values of Mach number M wa . Pr = 0.715; 7 = 1.400} 3 = 1) 
T vb* = !• 
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(a) log -J-p * against k. Dashed curves neglect terms. 

“ 0 d 

Figure 4.- Variations of zero-order pressure ratio with diameter ratio k 
for various values of Mach number. Pr = 0.715; 7 = 1.^00; p = lj 
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- First-order slip velocity at inner cylinder against diameter 
tio k. Pr = 0.715} 7 = 1.^0} 3 = 1} T w b* = 1} a = 1. 
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k 


Figure 7-- First-order temperature jump at inner cylinder against diameter 
ratio k.- Pr = 0.715; 7 = 1.400; 3 = 1; T w b* = 1; a = 0.900. 
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k 

Figure 8.- First-order temperature jump at outer cylinder against diameter 
ratio k. Pr = 0.715; 7 = 1.400j 3=1; T wb * = 1; a = 0.900. 





- First-order correction to parameter 
k for various values of Mach number 


against diameter 
a . Pr = 0.715* 


7 = 1.400* 0 


* = 1} a = 0.900} o=l. 





112 


NACA TN 2895 



0 -2 .4 .6 .8 1.0 

k 

i 

Figure 10. - 'First-order correction to pressure ratio 1 p b * against 
diameter ratio k. Pr = C.715j 7 = 1.400; 3 =1; T wb * = lj 
a = 0.900; a = 1. 




Figure 11.- Zero-order friction coefficient at inner and outer cylinders 
against diameter ratio k. Pr = 0.715; 7 = 1.400; 3 = 1; T wb * = 1; 

a = 0.900; a = 1. 




Figure 12,- First-order correction to friction coefficient for both 
inner and outer cylinders for various values of diameter ratio k 
Pr = 0 . 715 ; 7 = 1.400; p = 1; T w b* = !}<* = 0 . 900 ; a = 1. 
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k 


Figure 13*- Zero-order heat transfer at inner cylinder against diameter 
ratio k. Pr = 0.?15} 7 = l.^OOj 0 = 1> T wb * =. lj a = 0.900) a = 1. 
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k 


Figure 1^.- Zero-order heat transfer at outer cylinder against diameter 
.ratio k. Pr = 0.715; 7 = lAOOj 3 = 1; T w ^* = lj a = 0.900} a = 1. 




Figure 15 .- Variation of function 0.6336(1 - k) 2 (2 + k 2 )/2k 2 log ( 

against diameter ratio k. 





Figure 16.- Parameter r) against effective Knudsen number. 
7 = 1.400 j 3 = lj T w b* = 1} a = 0 . 900 ; a = lj k = 0.5; N 







cl> 
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Figure 18.- Slip velocity against effective Knudsen number. Pr = 0.715} 
7 = 1.1*00; P = 1} T wb * = 1} a = 0.900} a = 1; k = 0.5; M wa = 2. 
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Figure 19.- Heat transfer against effective Knudsen number. Pr = 0.715) 
7 = 1.1+00) p = 1) Twb* = 1) a = 0.900) a = 1) k = 0.5) M wa = 2. 
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Figure 20.- Temperature jump against effective Knudsen number. Pr = 0.715} 
7 = 1.400} 3 = 1} T wb * = 1; a = 0.900; a = 1} k = 0.5} = 2. ' 
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